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1. 


A non-Newtonian liquid is a liquid which does 
not conform to Newton’s postulate, which 
states that in simple shear the shearing (or tan- 
gential) stress p,; is proportional to the rate of 
shear (or tangential component of the velocity- 
deformation) é,:, on any plane whose normal is 
o wat = 

€:= ¢pr, 


(1)? 


where ¢ is a mechanical material-constant called 
the coefficient of fluidity, and the dot is used to 
denote the time derivative. This Eq. (1) is called 
the rheological equation of state of the material. 
It is a basic law and can as such not be directly 
verified.? If, however, the differential Eq. (1) is 
integrated for the stationary flow through a tube 
with the radius R and the length /7, the result is 


Q/t= eR*rp/8l, (2) 


where Q is the quantity of liquid flowing through 
the tube during the time ¢ under the pressure 
drop p/l. Eq. (2) is known as the Hagen’- 
Poiseuille law. If R and / are measured, and cor- 
related values of Q/t and p are observed, ¢ can 
be calculated from 


o=801/R*xtp (3) 


and is a constant‘ in the case of the Newtonian 
liquid if there is no ‘‘wall effect.’> For a non- 


1 The notation follows the system adopted in the author’s 
Introduction to Mathematical Rheology and explained in 
J. Rheol. 3, 264 (1932). Whenever there is no danger of 
ambiguity, we will, however, write p; and e for pr: and ént. 

2M. Reiner, Introduction to Mathematical Rheology. 
(In press); Cf. J. Rheol. 3, 245 (1932). 

3E. Hagen, Pogg. Ann. 46, 437 (1839) derived this law 
before Poiseuille. 

‘ With respect to all quantities appearing in Eq. (3). 

5 With respect to the “wall-effect” see Sections 8 to 11. 


TABLE I. Apparent fluidity of a solution of crepe-rubber in 
toluene (1.77 percent conc. by weight). Temperature 20°C. 








p* t 
cm H.O sec. 


10°Q/t ¢’ 


cm sec, rhes. 





(a) Measured in a Bingham viscometer 
R=0.0238 cm 1=9.93cm Q=5.98 cm? 











— 86.4 3038.1 1.31 1.23 
+131.5 1610.4 2.47 1.52 
— 181.0 1033.5 3.87 1.72 
+224.6 728.8 5.47 1.97 
— 274.2 545.0 7.32 2.16 
+316.8 425.0 9.38 2.39 
— 361.9 352.4 11.3 2.53 
(b) Measured in a Bingham viscometer 
R=0.0412 cm 1=9.90cm Q=4.00 cm? 

+ 38.3 846.8 4.73 1.10 
— 84.4 278.4 14.36 1.52 
— 132.3 142.4 28.1 1.90 
— 178.0 89.8 44.6 2.24 
— 225.6 62.4 64.2 2.53 








* The pressure entered here is the reduced pressure after taking 
Seegng cca neceneney canvectone. 2 pegiitre ppeeeweteagglind se te 
Newtonian liquid, however, which does not 
conform to Eq. (1) the quantity calculated 
according to (3) is not a constant and is then 
usually called the “apparent fluidity” yy’. (See 
Table I.) 

This statement cannot be reversed; i.e., if a 
material is pressed through a tube and a coef- 
ficient calculated according to (3) turns out to 
be variable, this does not prove that the material 
is a non-Newtonian liquid. It has first been shown 
by Bingham® that a material for which the 
postulate 


é.= o(p:— 8) (4) 


holds good, and which is a plastic solid with the 
yield value @, also leads to a variable apparent 


*E. C. Bingham, Bur. Stand. Sci. Paps. 13, 309 (1916). 
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fluidity, which is according to Buckingham’ and 
Reiner® 


= of 1—(4/3)(261/Rp)+(1/3)(261/Rp)*]. (5) 


For the non-Newtonian liquid therefore the 
criterion that gy’ is variable has to be amplified 
by the condition that @ vanishes. 


2. 


It is not possible to decide by appearance 
whether a given material which flows is a liquid 
or a very soft plastic solid. Every flow investiga- 
tion has therefore to be expressed at the first 
stage in terms of consistency which concept can 
be applied to liquids and solids as well. In every 
rheological experiment an observed kinematical 
quantity (V) is related to an observed dynamical 
quantity (P). In the case of the Hagen-Poiseuille- 
experiment the kinematical quantity is Q/t and 
the dynamical quantity is p or p/l. In the case 
of the Margules’-Couette” experiment, V is the 
number of revolutions of the outer cylinder per 
unit time (7) and P the difference of the readings 
on the deflection scale (As). In both cases however 
other more convenient quantities can be sub- 
stituted for the quantities just mentioned (see 
Section 6). If we plot the kinematical quantity 
as ordinate against the dynamical quantity as 
abscissa, the resultant curve is the ‘‘consistency 
curve.”’ (See Fig. 1.) 

From an examination of the consistency curve 
it is sometimes possible to decide whether the 
material, under investigation, is a solid. If the 
consistency curve starts at a point on the P-axis, 
other than the origin, the material is a solid. 
This statement cannot be reversed. If the con- 
sistency curve passes through the origin, the 
material may be either a liquid or a solid which 
does not adhere to the wall of the apparatus, i.e., 
which shows slippage or generally a wall-effect. 

It therefrom results that a theory of non- 
Newtonian liquids has to start from a general 
theory of flow of materials. At a later stage only 
(from Section 12 on) is it possible to branch off 
with a theory which is solely devoted to the flow 
of liquids. 

7 E. Buckingham, Proc. Am. Soc. Test. Mat. 21, 1 (1921). 

* M. Reiner, Koll. Zeits. 39, 80 (1926). 

*M. Margules, Wien. Sitzungsber (2A) 83, 588; 84, 49 


(1881) proposed the rotation-viscometer before Couette. 
10 M. Couette, Ann. de chim. 21, 433 (1890). 
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Fic. 1. Consistency curves of rubber-benzene solutions, 
the data for which are given in Table I ( plotted as 
usual; — - — - — plotted in consistency variables). 


3. 


A very valuable contribution to the general 
theory of flow has recently been made by 
Hohenemser and Prager.’ Every rheological 
equation of state is a relation between the tensors 
of stress and strain and their time derivatives. 
We confine ourselves here to homogeneous and 
isotropic or quasi-isotropic materials. We assume 
that the volumetric strain’? depends only upon 
the hydrostatic pressure and that vice versa a 
hydrostatic pressure does not change the form 
of the body, i.e., does not produce deformation. 
For the case of deformation, we therefore have 
to consider such stresses only as are in excess of 
the hydrostatic pressure and such strains as are 
in excess of the volumetric strain. It can be 
shown by mathematical analysis that in every 
such deformed body, where the volume has not 
been changed, there exist at every point direc- 
tions (n) for which the normal components of 
deformation (e,,) vanish, i.e., where there is 
only pure shear in the plane-element, which is 
at right angles to n. Likewise there exist in every 
point plane-elements the stress across which is 
reduced to a pure shearing stress. For an isotropic 
body it can be assumed that the plane elements 
of pure shear and pure shearing stress are iden- 
tical. With these assumptions, Hohenemser and 


1! K, Hohenemser and W. Prager, Zeits. f. angew. Math. 
u. Mech. 11 (1931); cf. J. Rheol. 3, 16 (1932). 

12 This is a strain which changes the volume only, but not 
the form of the body. 
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Prager have shown that if a relation between the 
tangential components of the stress and deforma- 
tion and their time derivatives is assumed, the 
rheological equation of state of the body is com- 
pletely determined; and they have also shown 
how the general three-dimensional differential 
equations, or the tensor relations, can be derived 
from the one-dimensional postulate about the 
tangential components. As a special case they 
have discussed the most general linear law of 
flow.'® It is however obvious that no linear law 
of whatever form can adequately describe the 
flow-properties of a non-Newtonian liquid. A 
theory of non-Newtonian liquids can therefore 
only result from a nonlinear law of flow. Such a 
nonlinear law will be derived in the next section, 
with Hohenemser and Prager’s system as a 
starting point. 
4. 


As has been explained in the previous section, 
the general rheological equation of state can be 
expressed by the formal equation 


F (bu Pu er &:)=0. (6) 


We are here concerned with a flowing material. 
In this case, the quantity e; cannot appear in the 
equation of state for the following reason: Flow 
is continuous deformation under a constant force. 
If the amount of shear is of influence on the 
state of stress in the material, this influence can 
only consist in increasing p, with increasing 
On the other hand, deformation under the 
action of external forces can only continue until 
the internal stresses have increased to such an 
extent as to counterbalance the external forces. 
When this state has been reached, the deforma- 
tion stops and consequently there is no flow. In 
a flowing material, which may be either a liquid 
or a plastic solid above the yield value, the 
amount of deformation is therefore not material 
but only the deformation velocity. This is tan- 
tamount to the stipulation that we do not con- 
sider strain-hardening. 

For a flowing material Eq. (6) becomes accord- 
ingly: 


f(Pu pu &)=0. (7). 


Eq. (7) can be further simplified if we consider 


13K. Hohenemser and W. Prager, J. Rheol. 3, 16 (1932). 
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steady flow only. This is the case with the 
Hagen-Poiseuille, and the Margules-Couette 
experiments which we may call the fundamental 
rheological experiments. If the flow is steady, 
the stresses are for obvious reasons functions of 
the space coordinates only, but not of the time. 
Therefore p, vanishes and for stationary flow we 
have 


f(t, €:)=9. (8) 


This result is of importance for the following 
reasons: In Europe ‘‘viscosity anomalies” have 
first been observed in the case of liquids showing 
dynamic elasticity, i.e., elastic behavior against 
instantaneous stress.'4 This leads to explaining 
these anomalies as due to an “‘elasticity”’ of the 
material. This elasticity was then identified with 
the static elasticity of cohesional forces and a 
law of the same form as Bingham’s law (4) 
assumed by Szegvari® and Reiner.’ The law, 
Eq. (4), however, covers a wider range, as no 
assumption is introduced with regard to the 
origin of the yield value @. The yield value may 
be due to cohesional forces showing elasticity; 
it may, however, also be due to other forces, as 
solid friction between suspended particles, sur- 
face tension, etc. As a matter of fact, Bingham 
had first, when proposing his equation, started 
from a concept of solid friction and called 6 
accordingly ‘‘the friction constant.’ It should 
now be remembered that the phenomenon of 
dynamic elasticity is connected with the oc- 
currence of p, in the differential equations of 
flow. As it has just been proved that in the 
case of stationary flow ; vanishes, this also 
proves that there is no direct connection between 
an observed dynamical elasticity of a material 
and its viscosity “anomaly” in either of the 
fundamental rheological experiments. 


5. 


Instead of Eq. (8), we can also write 
é:=f (pr). (9)" 


The function f can be developed into a Taylor- 


14 W. R. Hess, Kolloid Zeits. 27, 154 (1920) where refer- 
ence is made to observations of 1906. 

18 A. Szegvari, Zeits. f. physik. Chemie 108, 175 (1924). 

16 M. Reiner, J. Rheol. 3, 253 (1932). 

17 Here f indicates, of course, a different relationship than 
f in Eq. (8). 
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series and the result is 
amo 1 df(™ (9) 
hit TE an eenigen DP. (10) 
n=in! dp," 


In Eq. (10) we have an infinite number of coefficients {‘"(@). They are, however, all functions of a 
limited number of physical constants, which characterize the rheological behavior of the material. 
Eq. (10) is the general law of steady flow of any material. It was first derived by Reiner" without, 
however, noting the restriction that it holds good in the case of steady flow only. 

The problem is now, given the general law of steady flow, Eq. (10), to find the corresponding V/P 
relation. The procedure is quite analogous to the one adopted by the author in a previous paper." 
By introducing in the case of the 


Hagen-Poiseuille Experiment Margules-Couette Experiment 
pr=rp/2l, p:=(}xr*)(M/I), (11) 
é,= —dv/dr, @,=rdw/dr, (12) 
where where 
p/l is the pressure drop, M is the rotational moment, 
1 is the length of the capillary, l is the length of the inner cylinder, 


r is the distance of the layer from the axis of the cylinder, 














we get 
dy _*y £2) (1? )’ # 2p LO) ( M _ y 
~ dr po (55 ee -_ a 2rrl o)» (13) 
and through integration 
_ 2b" £0) (rp/2l—o9)yn"" . £9)" | "FO (- 1)" ( n ) (=e) 
ar myil a ae ~, Um—n)\m)\ 
for pr=d, (-1)"1 k 14) 
29 +(-1)"Inrp +k ( 
i.e., for r=ro= = 
for pr=d, 
i.e., for r=ro=(M/2rid)!, (15) 


For the determination of the integration constant k it is usually assumed that the velocity of the 
flowing material in the layer adjacent to a solid wall is identical with the velocity of the wall. If 
this is the case we have 





pa — 2S £0) (Rp/21-0)"™ pa 3S Lo” mm 1)" (” )( M y" 
nm os n! n+1 z aan n! foe 2(m—n) \ m 20K) 
where R is the radius of the capillary, +(—1)" In Ri ; (16) 


where R, is the radius of the inner suspended cylinder, 


and therefore 











2S (8) (Rp/21—3)"*! —(rp/21—9)"™" on T LOW)” "3 {< (" ) ( zh i 
= hed n! n+1 7 =“ n! ~ a(n me) \ aad 
for r=So, 


x[R? m—n) — 2(m—n)] + ( —])n In c/o} (17) 
for ro, 


18 M, Reiner, * J. Rheol. 1, 11 (1929); > 1, 250 (1930); °1, 5 (1929). 
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Fic, 2. Corresponding flow conditions in capillary tube and rotation instruments. 


while for 


r—TPo 





f—ve 
we have 
pave 2E"E £0) (Ro/2 dn | sone 2 Oe "Ss" 
- p bea n! n+1 . jer! n! | ee 2(n—m) 


M \"™ M_\ 
xX [(<m) -1| +(—1) in( =a) . (18) 


See Fig. 2. 


The volume of flow is therefore 


v7=R 


O/t=rétxvot2x f vrdr 
T=To 


_ 82/3 oy f™(9) 


PX, (n+3)! 





(Rp/21—a)"*" 


< [(Rp/21)?(n +1)(n+2) +(Rpd /l)(n +1) +287). 


The angular velocity of the outer cylinder with the radius 
R2 is therefore 


ee Fem(ayan "et (1) (" ) ( M 


1 (n— 
“ae n! | mao 26—m) \m 2nld 





X[R2e"-™ — R,2"-™) }+-( —1)" In R/R,| (19) 


for M=2rldR:*; 
and 


n=O | man—1 


- f£™(9)0" = (—1)* [( M " -1| 
eo! n! | aa 2(n—m) 2rldR;? 


mM \}| 
—_ n nmin: 9g’ 
+(-—1) in ( 5533) (19’) 
for M=2rldvR?. 


Q= 








Eqs. (19) express the general V/P or consistency relations for both fundamental rheological experi- 
ments. This solves the problem which we put before ourselves in the beginning of this section. If the 


Eq. (10) is given, we know all the coefficients f’(#) 


,f'(d) ++ +f (9). We therefore also know the coef- 


ficients of the power series Eqs. (19) and (19’). Conversely, if the consistency relations are empirically 
given through observation and the consistency curves can be analyzed, so that they can be ex- 
pressed through power series of known coefficients, this leads to the establishment of the general 
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Fic. 3, Consistency curves for three types of material. 


law of steady flow, Eq. (10), for the material under investigation. As will be shown later (Section 14) 
this method is however very limited in its possible applications. 

If we specialize Eqs. (19) and (19’) for n=0, we get the equations for the non-Newtonian liquid; 
if we specialize for n=1 we get the equations for the Bingham body; if we specialize for 8=0 and 
n=0 we get the Hagen-Poiseuille and the Margules-Couette laws, respectively. 


6. 


In the form just given, the consistency curves are not independent of the dimensions of the ap- 
paratus. (Compare Fig. 1.) It should however be noted that neither is this the case with, e.g., the 
Q/t : p/l relation of the Hagen-Poiseuille law. General considerations of similitude lead us to expect 
that it should be possible to select such expressions for the V/ P variables that the consistency curves 
for the same material and similar apparatus will be independent of the dimensions of the apparatus 
irrespective of the laws of flow and also of whether the material is a liquid or a solid. 

An examination of Eqs. (19) together with expressions (11) shows that if we plot 


V’ =40/rR°t V=20/(1—a), 
(20) 
a=(R,/R:)? 
against the shearing stress at the wall 
P’=Rp/2l P’ = M/2rlR,;?, (21) 


we have in 


, 1 —— 4f(™ (9) 7 
V’= > ofa  _ a\n+1 2 as - 
Pon (nay) P98)" +1)(n +2) 


1 "3 fo~r(ayan { "1 (1) ( ) 
i-a 2! | nap 2am \m 
+2(n+1)3/P’+2(8/P’)*] ; 
am Pas x(P/d)"-™[1 —a"-™]—(—1)" Ina} (22) 


for P’'=d/a 
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and 
Fr ie /— ; i 
J =0 for I =v, a 1 "Sf (a) 9" . "(1)" ‘ ) 
l-a on! np Bam \m 
x [(P’/8)"-™" —1]+(—1)" In P’, o (22’) 
for P'=3/a 
and V’=0 for P’=d, 


expressions which are independent of the dimensions of geometrically similar apparatus." (See Fig. 1.) 
Such a demonstration as given here for the fundamental rheological experiments is always possible 
and no other V’/P variables should herefore be used in plotting consistency curves. For this reason, 
the concept ‘‘consistency curve” will be used so as to mean the curve when plotted in such a way as 
to be independent of the dimensions of the apparatus and we will write V and P for V’ and P’. (See 
Fig. 3.) 
va 
The apparent fluidity can be calculated from expressions (22). We get 





ga 80 VL i a ti a ,Z + 
Rip PO Pt M(1/R2—1/R2) (1—a)P 
for P=8 for P=3/a 
and g= ant y= I Pp = «ee (23’) 
for PSs. (1—a) n=l 
for P=3/a 
and ¢y’=0 for P=d. (23’’) 


If we plot the apparent fluidity against the shearing stress at the wall, the resultant curves are 
therefore also independent of the dimensions of the apparatus. 
If we now examine the equations of state (1) for the Newtonian liquid ond (4) for the Bingham- 


body, it can be seen that the quotient é,/p, gives the fluidity ¢ in the first case, but not in the second 
case. In the second case we get 


é; ‘pbi= o(1—0/p,), (24) 


and this is the apparent fluidity for the plane two-dimensional case. (See Fig. 4.) The fluidity can 
therefore generally not be defined by the quotient of the rate of shear to shearing stress and different 
definitions have to be used in the case of the fluid and the solid for a magnitude which is of the same 
nature. For this reason also a different name, ‘“‘mobility,”” was introduced by Bingham in the case 
of a solid for what is “‘fluidity”’ in the case of a liquid. 

If, however, we define ¢ by 


g=dé,/dp. (25) 


we get in both cases, i.e., from either Eg. (1) or (4) dé,/dp,=y. We are then justified in using one 
term only, i.e., the term fluidity. In analogy to the expressions (25) we define by 


=dV/dP (26) 


the practical fluidity ¢*, where V and P are such kinematical and dynamical quantities as make the 
consistency curve independent of the dimensions of the apparatus. 

The usefulness of this new concept will be seen from the following: Suppose we are making flow 
experiments with a material of which we do not know whether it is a liquid or a plastic solid, we 





1° For the tube this qualification is not necessary, for the 20 Compare with the usual nomenclature, 


e.g., Mooney, 
rotation instrument however this implies that aisconstant. J. Rheol. 2, 215 (1931). (Editor.) 
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Fic. 4. True and apparent fluidity of the Bingham body. 


measure related values for 
Q/t and p/l 2 and M 
calculate - 
"=40/rR% and P=Rp/2l V=20/(1-—a) and P=M/2rR,’l 


from the known dimensions of the apparatus”! and plot V against P. From this curve, we determine 
graphically (or by some other suitable method) ¢* and plot it against P. This curve then represents 
the functional relationships: 








— —_ 4f(™ (3) - ee mm Fin) (9) oe - _— 
e" P bol (nai?) (n(n+1)(n-+2) en fe bod n'!P [(P—9)"—(Pa—9)"] (27) 
+3n(n+1)3/P +6n(3/P)?+6(8/P)*} an P=8/a 
for P= neco 
a —.. —— a 
and ¢*=0 and ¢ T=a 2 nIP (P—0) (27’) 
” “ae for P=38/a, 
and ¢*=0 (27”) 
for Pd. 


For the Newtonian liquid, i.e., 2 =1, 3 =0, specialization of the Eq. (27) gives in both cases 


¢* =f'(0) =¢. 
For the Bingham-body, i.e., »=1, we get 


¢* =f'(9)(1 —(8/P)*] = ¢[1 —(8/P)*] g*=f'(9)=¢ for PZb/a (28) 
> ‘(9 ’ 
for PE), and t= 2) (1-9/P)=[e/(1—a)L(1-9)/P] (28) 
which for large values of P can be put as equal to ¢, ” _ 
for P=8/a. 
*1 Instead of determining directly the dimensions of the 2 After differentiation of the expressions on the right 


apparatus, we can of course calibrate it with a liquid of | side, Eqs. (22) and (22’), the m-summation can be carried 
known constants, out. 
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Fic. 5. True, apparent, practical, and mean practical fluidity of the Bingham body. 


while for P= 3, we have ¢* =0 in both cases. The practical fluidity therefore has a very close relation 
to the true fluidity; as a matter of fact, both are for all practical purposes identical. In Fig. 5, the 
true fluidity, the practical fluidity and the apparent fluidity are shown and the advantage of using 
the practical fluidity is obvious. Besides the true and the practical fluidity, the following quantity 
is of interest: 


m= é:/(pr— 9). (29) 


We may call it the mean fluidity. This designation suggests itself from an inspection of Fig. 5. The 
mean fluidity ¢,, for the point (p;, é,) is identical with the true fluidity ¢ of a point between (0, 0) 
and (f;, é,). In the case of a liquid 3 vanishes and the apparent fluidity ¢’ is identical with ¢,,. We 
will accordingly use the following terminology. The apparent fluidity 


= rate of shear y « kinematical consistency variable 
Fee , RCMarical CONE 3 


y= e. or i.e. _—__—___——_——— 
pe shearing stress Fr’ ’ dynamical consistency variable 


The apparent fluidity will only be used if it is not known whether the material is a solid or a liquid, 
or in case it is a liquid showing a ‘“‘wall-effect.” 
True fluidity p=dé,/dp,, 
mean fluidity yg, = é,/(p:—#), 
practical fluidity g*=dV/dP, 
mean practical fluidity y,,*= V/(P—P»), 
where Pp is the dynamical quantity corresponding to the yield value or the abscissa of the “yield 
point.” 
Besides these quantities it is sometimes advantageous to use another quantity, which, without 


being a “‘fluidity” is very closely related to it. Reynolds has introduced the plotting of flow-curves 
on a logarithmic scale. If we put 


V.=log V, P,=\log P (30) 
the V; : P; curve is the logarithmic homologue of the consistency curve. We can then define by 
¢gi=dV,/dP, (31) 


the logarithmic fluidity. The usefulness of this concept will be shown later. 
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8. 


From Eq. (27) we find that generally ¢* vanishes for P = 3. This could therefore furnish a criterion 
whether a given material is a liquid or a solid. As has, however, been remarked in Section 2 this 
criterion holds good only where there is no ‘“‘slippage”’ or generally no ‘‘wall-effect.”’ 

As Bingham has pointed out, it is ‘quite possible that the fluidity of a liquid near the surface 
(solid wall) is not identical with that within the body of the liquid.”’ ** This gives rise to a wall-effect. 
However, even if the material does not wet the solid wall, one cannot properly speak of slippage in 
the sense in which this word is ordinarily used. In this case, there exists, according to Bingham,”* a 
layer of some other medium, usually air, between the material and the solid wall. In the case of a soft 
solid, which consists of a liquid and a solid phase, as e.g., an oil-paint consists of the liquid vehicle 
and the solid pigment, a wall layer of pure liquid has to be assumed to be present between the 
wall of the apparatus and the flowing solid (e.g., the oil-paint). 

We have therefore to generalize our equations to include such cases. Let us denote the thickness 
of the wall layer by d and its fluidity by ¢,. We assume ¢,, to be constant, i.e., the wall layer is a 
Newtonian liquid. 

The wall layer adheres to the wall.2* We therefore have 


V = ¢u(p/4l)(R?—r*) w= go( M/2nl)(1/R2—1/r) (32) 
for the layer near the wall. for the layer near the inner cylinder, and 
w= gy(M/2rl)(1/r? —1/R:2?) (32’) 


for the layer near the outer cylinder. 


These equations follow from Eqs. (17) for #=0, »=1. The velocity at the interface of the w non layer 
with the body of the liquid is 





d = gu(p/4l)(R? —(R—d)?*]. ‘iets a ( J 7 —? (33) 
M 1 1 ' 
04," Pv 5. ni (ca ~ Re)” (33’) 


(See Fig. 6.) We now have to go back to Eqs. (14) and to determine the integration constant so that 





gat "sy £20) [(R—d)p/2l-—3)}""" 7 ‘ _"s £(o)o" — (—1)" (”) 
ions n! n+1 — “1 en n! 9 2(m—n) \m 


M —m : 
x (sg says) +(—1)" In (Rt) +h (34) 


Introducing these integration constants into Eqs. (14) we get the following equations instead of 
Eqs. (17) 














_2l "> S™ (8) [(R—d)p/2l—8 }"*' —[rp/2l—8]}"* _ = f™(9)d" rs * (-1)" (")( M ~ 
p pool n! n+1 sagt d n! | inal 2(n—m) \m/ \2nld 
D «ms P 
+¢e 5 [R'—(R-d)*] Lasdie~siabit- 
4] X| (Ritd an +(—1)" in 75 
for (R—d)=r=ry=201/p. 


Tew aa Ge Trap) (35) 


for (Ri +d) SrSro=(M/2ai0)!. 
Furthermore there is 
Vw = ¢w(p/4l)(R? —r*) Ww = Gw(M/2rl)(1/R2—1/r) (35’) 
for R=r=(R—-d) for RiySr=(Ri+d) 


*%E. C. Bingham, Fluidity and Plasticity. * M. Reiner, Zeits. f. Physik 79, 139 (1932). 
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2e.P | 1— : -| (36 
+200? [\-GaaR sa RF) 6°) 
for P=8/a, 
and 
y 1—2d/(R2—R)) ] ” 
=29,P}1— - - 36 
es [: (1+d/Ri)*(1 —d/R2)* sii 
for P=8(1+d/R;)?. 


From Egs. (36), (36’) and (36’’) it follows that whenever there exists a ‘‘wall-effect,” the con- 
sistency-curve is not independent of the dimensions of the apparatus and furthermore that, if the ma- 
terial is a plastic solid, the consistency curve starts with a straight lirfe. The second part of this 
statement cannot be reversed, i.e., if the consistency curve starts with a straight line, the material 
may either be a plastic solid with a wall-effect or a liquid with or without a wall-effect. 


9. 


Expressions for the fluidities in case where there exist wall-effects can be easily derived from Eqs. 
(36) and (36’). For the Newtonian liquid, i.e., n=1, 8 =0, specialization of the resultant expressions 


gives 


1—2d/(R2—R 
¢* = gu—(¢u—¢)(1—d/R)4. gt =2 /(Ro 1) 


+d/Ry¥A=d/Rp = 8) 





¢w—(2¢~—¢) ai 
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Fic. 6. Plastic flow in the case where there exist wall-layers. 


For the Bingham-body, i.e., n=1, we get 
¢* = gw—(¢u—¢)(1 —d/R)*— (0 P)* 1—2d/(R2—R:) (38) 











¢* =2eu—(2¢eu—¢) 1/R 2(4 d R 2 
for P=#/(1—d/R), “ (1+d/Ri(l —d/R2) 
for P=I/a 
and o* = ¢u[1—(1—d/R)*] and 
for 3/(1—d/R). a [1- 1—2d/(R:—R)) | 
a (1+d/R,\)2(1 —d/R2)* 
1 1 v , 
+4 ele (38') 
for P=38/a, 
and 
1 —2d/(R2—R) | = 
p* =2eul 1— 3 
” e-[1 (1+d/R,)*?(1—d/R2)* om 
for P=o0(1+d/R,)?. 


For the non-Newtonian liquid, i.e., ) =0, we get 











eee ee pated 1 —2d/(R>—R;) 
et age—toe— 5 LO PU-d RP"! y_ayry. = 2%! — Tak —d/ RP 
[ ome n+3 (n—1)! : 
= For 1 et | (39) 
‘l-e n! (1+d/R,)" (1—d/R2)" ; 
; 10. 


It was stated in Section 2 that if the consistency curve of a material starts from a point on the 
P-axis with a finite abscissa, this proves that the material is a solid. We are now in the position to 
amplify this criterion for the case where the consistency-curve passes through the origin. In this 
case one has first to determine the practical fluidity at rest, i.e., for P =0. If the material is a plastic 
solid this is 





1—2d/(R2—R) | 
°=g.1—(1— ‘ eT ; 
¢o0* = gw 1—(1—d/R)*]. ¢o 2ee| 1 (i+d/R,y° —d/ Re) (40) 
If however it is a liquid we have 
¢o" = got¢uL:: ]. ¢o* = gotle¢uLl:::]. (41) 
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If we now vary the dimensions of the apparatus, i.e., 
the radius R the radii R,; and Ro, keeping a constant 


and plot the calculated values against 
1/R 1/R, 
then extrapolation of this curve for 


1/R=0 1/R,=0 


gives go* =0 in the case of a plastic solid and yo* = go¥0 in the case of a (Newtonian or non-New- 
tonian) liquid. 


This provides us with a definite criterion of whether a. material is a plastic solid or a liquid. As 
we are not, in the present paper, concerned with the flow of solids, we will assume in our further 
investigations that #=0, i.e., in the general case of a wall effect we have 


V=Pedi—(1—d/R)] ~_ ie.) 
"E 4620) P"(n+1)(n+2)(1—d/R)"*8 ¢ (OCs, mains 

=, (n+3)! . 1 "SS f~™ 0) = 1 _ a” | (42) 

+ . (i+d/Ry™~ (—d/Ry™ | 


l—a n! n 


n=l 





V=2P en [ 1 








and where there is no wall-effect we have 


— 4f(™(0) ™~ (0) Pr 


, ' : p 1 ee. | el 
V= pers (nt3yi? (n+1)(n+2). V=  — —a n! - (1 a ;™ (43) 


These consistency curves correspond to a rheological equation of state of a non-Newtonian liquid 





. nao f(") (Q) 
é:=f(pi)= ~ a! pr”. (44) 
11. 


The wall-effect can be considered under two aspects. It may be of interest to determine the thick- 
ness of the wal! layer and its fluidity. In order to accomplish this, we differentiate Eqs. (41) in 


respect of ’ 
1/R 1/R, remembering that 
R,/R2 is a constant. 
We get 
dgo*/d(1/R) =4¢d(1 —d/R)*. dgo*/d(1/Ri) =4¢ud 


((14+d/Ri)(1—d/R2) +R2/Ri-1 —2d/Ry_ 
(R2/Ri—1)(1+d/R,)*(1 —d/R2)§ 





(45) 
Extrapolating the curve defined by (45) for 1/R=0 we get in both cases 
Ldgo*/d(1/R) Jo= 4¢ud. (46) 


If ¢» is known this gives the thickness of the wall layer. Otherwise, a further differentiation and 
extrapolation furnishes another value so that both unknowns ¢, and d can be calculated. 

If, however, the problem is to determine the consistency of the body of the material, the wall- 
effect is bothersome and has to be eliminated by means of correction. This can be done as follows: 
One has first to determine the true fluidity of the material at rest go as explained in Section 10. If 
we know this, we redraw the consistency curve for a coordinate system as shown in Fig. 1. The 
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consistency curve in respect to the new VP coordinate system is given by the following equations: 
[ V— P(¢o* — ¢o) leu [ V— P(¢o* — ¢0) 2 eu 
Gut ¢o—¢o* 2¢wt+¢o— ¢o* 
"> 4f(0) sii "y £() Pr (1 —d/R2)*" —a(1+d/Ri)") 
= >» Pr» 2 onde} n 1, = >. < a 
2 Gta OTE +20 —d/R) zat nm ((—d/Ry? all +d/R7] 


1 
" 1+d7RiP" (1 —d7 Ry 

















(47) 


The first term of the series, i.e., for »=1, is f’(0)-P and is identical with the first term of a liquid 
showing no wall-effect. In the case of a Newtonian liquid showing a wall-effect this curve is therefore 
the correct consistency curve of the bulk of the material. We may therefore call the consistency 
curve as given by Eq. (45) the “‘bulk’’-consistency curve. It is identical with the consistency curve 
of the bulk of a non-Newtonian liquid showing a wall-effect, if either small P’s are considered or 
d/R is small. For larger P’s and larger d/R’s there is no way of plotting consistency curves inde- 
pendently of the dimensions of the apparatus, if the law of flow of the material is not known. 





12. 

The foregoing investigations have furnished 
us (a) with a criterion with which to decide 
from observations whether a material, the law 
of flow of which is not known, is a solid or a 
liquid, and (b) with a method of how to plot the 
observations so that the curve is entirely inde- 
pendent of the dimensions of the apparatus, 
where there is no wall-effect, or is approximately 
independent of the dimensions, where there is a 
wall-effect. 

It is, however, not possible to plot the con- 
sistency curve in such a way that observations 
made with the same material in different sort of 
apparatus (e.g., capillary and rotation) will give 
the same curve, if the law of flow of the material 
is not known. The only advantage, at this 
stage, of having made observations with a 
different sort of apparatus is this, that the 
different consistency curves have the same shape 
near the origin or all extrapolate to the same 
practical fluidity at rest, which is identical with 
the true fluidity at rest. It is therefore possible 
to compare the measurements made with dif- 
ferent sorts of apparatus near the origin without 
the law of the material being known. The next 
and more difficult problem is therefore to find 
the law of flow of the material. 

Two methods are possible. They have been 
named by M. D. Hersey* the differentiation 
method and the integration method. The first 
consists in the following: After having plotted 


25M. D. Hersey, J. Rheol. 3, 23 (1932). 





the consistency curve, let us assume that it is 
possible to express it in the form of a function 


V=F(P). (48) 


Then it has been shown by Weissenberg”™ in the 
case of the capillary instrument and by Mooney” 
in the case of the rotation instrument, how it is 
generally possible to derive from (48) through a 
differentiation process the corresponding Eq. (9) 
for the material under investigation. The dif- 
ficulty, however, lies in finding Eq. (48). Weis- 
senberg’s method has actually been applied only 
in the case of simple polynoms as 


V=a,P+aP? (49) 
or 


V=a,P+aP3 (49a) 


and these cases are special cases of the power 
series Eqs. (43). There exist different methods 
for determining the coefficients of the power 
series Eq. (43), and if these are known, the coef- 
ficients of the power series Eq. (44) are also 
known and therefore also the law of flow of the 
material. 

One method consists in extrapolating the 
V-curve and successive “‘slope’’-curves for P =0. 
Another method which is based on the plotting 
of the logarithmic homologue of the consistency 
curve has been used by Weissenberg, Rabino- 


26B. Rabinowitsch, Zeits. f. physik. Chemie A145, 1 
(1929). 

27 M, Mooney, J. Rheol. 2, 210 (1931). 

28 Such an equation was previously proposed by R. V. 
Williamson, reference 30. 
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witsch and Eisenschitz.?® These methods, which 
are essentially graphical methods, are found in 
actual practice to be too uncertain. The uncer- 
tainty of extrapolation is well known and the 
logarithmic method is in this respect not better. 

With a limited number of terms it is, of course, 
possible to calculate the coefficient a to any 
degree of accuracy by introducing the observed 
pairs of values V/P and solving the linear 
equations for the unknown coefficients. For an 
infinite series the application of this method is 
only possible if the series converges quickly and 
can therefore be approximated by a limited 
number of terms. As will however be shown 
later, the convergence of the series (43) is usually 
not very good and as to rheological equations of 
state of the form (49) or (49a), ie., with a 
limited number of terms, there exists a serious 
objection against them, with which we also will 
deal later. 

13. 


The integration method consists in postulating 
a law of flow of a hypothetical material and 
checking on the observations whether such a 
material is a proper model for the material under 
investigation. If such a law @,=f(p;) is assumed, 
it has to be integrated in order to obtain the law 
V=F(P). If this cannot be effected in closed 
form, f(p:) is best expanded into a power series 
and the integration then leads to a series of the 
form of Eqs. (43) [or (42) ]. This method has 
been applied successfully in the case of the 
Newtonian liquid, the Bingham body, etc. The 
differential equations which form the postulate 
for these model substances furnish when in- 
tegrated formulas which quite adequately de- 
scribe the flow properties as observed with large 
classes of material bodies. The laws suggested 
themselves by their simple form: they are linear 
laws. But no linear law is suitable to describe 
the flow properties of a non-Newtonian liquid. 
Following this clue of ‘‘mathematical simplicity” 
the problem would therefore be to assume a 
basic law of as simple a form as possible which 
conforms to the general shape of the consistency 
curves of non-Newtonian liquids, to integrate 


this basic law-differential equation, and compare 


29 R. Eisenschnitz, B. Rabinowitsch and K. Weissenberg, 


— Materialprufiingsanstalten 9, Sonderheft 91 
). 


the resulting consistency curve with actual 
measurements. 
14. 

If we examine the consistency curves of non- 
Newtonian liquids as observed and published 
until now, we can say that they generally show 
the following characteristics: The curve starts 
at the origin and the ordinates (V) increase 
monotonously with the abscissa (P), i.e., there 
exists no maximum JV. In addition to this the 
curve generally turns its convex side towards the 
P-axis. With regard to small V/P values many 
investigators found that (a) the curve has a 
finite slope at the origin. With regard to the 
shape of the curve for large V/P values, Ostwald 
and his school maintained that (b) there exists 
an inflection point, after which the curve gradu- 
ally becomes a straight line which, by extra- 
polation, passes through the origin (compare 
Fig. 7). Apart from these empirical findings one 
can say a priori that while the characteristic 
property of a non-Newtonian liquid is given by 
the fact that its fluidity increases with increasing 
shear, this fluidity naturally cannot increase 
beyond any limit and must finally reach a 
limiting value. We may denote this value by ¢,, 
bearing in mind that it may of course be reached 
at a large, but finite shear. 

In the light of the condition (c) it is clear that 


no polynom with a limited number of terms such 
as 


V=a,P+a2P*+-+-+a,P", n+ (50) 


can form the rheological equation of state of a 
non-Newtonian liquid. For if we divide V by P, 
we get 

gy’ =a,+a2P+:-++a,P"", (51) 


and unless 
de=a3= +++ =a,=0, (52) 


which is the case of a Newtonian liquid, 9’ 
cannot have a finite value for P= ©. This is the 
objection which rules out rheological equations 
of state of the form of Eqs. (49) and (50), or 
generally (51). y’ can, however, have a finite 
value for P=0 if m in Eq. (51) is infinite. The 
power series so formed, which is identical with 
the power series Eq. (23) for 6=0, is the develop- 
ment of a function which has a finite value for 
P=. In special cases, i.e., if the radius of 
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Fic. 7. General shape of the consistency curve of a non- 
Newtonian liquid. 


convergence of the power series is infinite it 
may even be possible to find the value for ¢,, 
from the series by summation without the func- 
tion f being known, and if the series converges 
quickly such an expression as Eq. (51) might be 
a suitable approximation of the series. In order, 
however, to be justified in taking only two terms, 
as has been done by Rabinowitsch, or three 
terms, as has been attempted by Reiner, the 
series must converge very quickly. An examina- 
tion of the papers mentioned under reference 18 
shows that in these cases this was not the case. 
A polynom can therefore not form the rheo- 
logical equation of state of a non-Newtonian 
liquid. Only in a limited region in the vicinity of 
the origin can it form an approximation ex- 
pression of the equation of state. Generally a 
polynom is therefore not suitable. Looking for 
another simple mathematical expression, it can 
be shown that there exists an algebraical function 
which conforms to all the conditions enumerated 
at the beginning of the previous section. Such 
a function has been proposed by Kraemer and 
Williamson™ and by Reiner.* The objections are 
however that its coefficients do not allow of a 
physical interpretation and moreover when in- 
tegrated the algebraic function becomes a tran- 
scendental one and thereby loses its only 
recommendation, i.e., simplicity. 


30 E, O. Kraemer and R. V. Williamson, J. Rheol. 1, 76 
(1929). 
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15. 


We are therefore forced to abandon the idea 
of a simple algebraical expression. This, however, 
robs us of our leading principle in our search for 
the law of flow of the non-Newtonian liquid. We 
are accordingly forced to look after some other 
indication of how to proceed. We may hope to 
find the right way in which to proceed if we 
abandon the purely formal treatment and form 
for ourselves physical concepts through which 
the consistency curve of a non-Newtonian liquid 
may be interpreted. Incidentally, we have 
already become acquainted with two such 
concepts, i.e., go the fluidity of material at rest 
and ¢,, the maximum fluidity which the material 
approaches at infinite shear.*! It should be noted 
that these quantities are the same whether we 
consider the apparent or the practical mean or 
true fluidity. The reason is that all liquids are 
comparable at the states of rest and infinite 
shear, i.e., P=0 and P= are comparable states. 

We can now say that the fluidity of the non- 
Newtonian liquid increases with increasing shear 
from ¢» to ¢,,. The flow condition of the non- 
Newtonian liquid at a certain shear is therefore 
characterized by the value of »9,,—¢. For a 
Newtonian liquid this value is equal to zero. The 
increase from ¢ to ¢,, or the consummation of 
n—g can be more or less rapid, according to 
the value of dy/dp;. Between both these quan- 
tities there exists the relation 


C=(¢,,—¢)/de/dpr. (53) 


C is the measure of the hydrodynamic “‘stability”’ 
of the liquid. If C=0 (see Fig. 8) the viscosity mo 
of the liquid at rest is at once reduced to 7,, as 
soon as the liquid is sheared. The liquid then 
behaves as a Newtonian liquid with the constant 
viscosity 7,, or the constant fluidity ¢,. If C= 
the viscosity mo of the liquid at rest is not changed 
at all through shear and in this case also we have 
a Newtonian liquid with the constant viscosity 
no or the constant fluidity go. We may therefore 
call C the ‘“‘stability coefficient’ of the liquid 
or 1/C its “lability coefficient.” 

A special class of liquids would be the class of 
all liquids for which C is a constant. If this 
constant is finite, the liquid is a non-Newtonian 
liquid. 


31 Care is taken to avoid turbulence. 
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If C is a constant, Eq. (53) can be integrated. 
The result is 


P= Poo — (Peo — Poe Pe! ©. (54) 


An equation analogous to Eq. (54) can be derived 
by considering a function relationship of (¢«) 
instead of ¢(p,). The result is 


n= Not (no— nn)e~*t! ©. (55) 





and 





For p,=0 and p:=® both equations give the 
same results ¢,,=1/n,, and go=1/m because in 
both these cases @; is also 0 and ©, respectively. 
Generally however ¢ as calculated from (54) is 
not equal to 1/7 as calculated from (55) and we 
have two different laws. 

It is possible to integrate Eq. (54) in closed 
form. We get 


V= 9,P —(24C*/P*)(¢.— 90) {1—e-?!° + [14+ P/C+ (1/2)(P/C)?+ (1/6)(P/C)* J} (56) 


g*= 9+ (72C*/P*)(¢.— go) {1—e-? [14+ P/C+(1/2)(P/C)?+ (1/6) (P/C)§+ (1/18)(P/C)* J}. (57) 


Differential Eq. (55) cannot be integrated in 
closed form. An attempt at an integration 
through development into a series was made by 
Reiner and Riwlin.® Eqs. (56) and the analogous 
equations derived from Eq. (55) give consistency 
curves which conform to the general shape of the 
consistency curve of a non-Newtonian liquid as 
described in Section 14. 


16. 

Eqs. (56) and (57) enable us to correlate 
measurements made with different sorts of 
apparatus as explained in Section 12. Moreover, 
the coefficients appearing in the rheological 
equation of state allow of a physical interpreta- 
tion in terms of macromechanical quantities. 
This is as far as the phenomenological point of 
view can lead us. It does, however, not enable us 
to gain an insight into the processes in the interior 
of the liquid, which are the cause of its non- 
Newtonian behavior. In order to gain such an 
insight, we have to proceed from the phenomeno- 
logical point of view to the structural, by con- 
sidering not only macromechanical, but also 
micromechanical processes. 

The first thing then is to note that all non- 
Newtonian liquids are dispersion systems. No 
pure homogeneous liquid is known which would 
not follow Newton’s law. The non-Newtonian 
behavior of a liquid must therefore be due to the 
changes which are produced in the dispersed 
phase through the flow. These changes may be 
either changes of the primary particles them- 
selves or changes in the structures of secondary 


32M. Reiner and R. Riwlin, Koll. Zeits. 43, 72 (1927); 
ibid. 44, 9 (1928). 





particles formed by the primary particles. The 
view that the flow anomalies of the non-New- 
tonian liquids were due to structural changes was 
first suggested by E. C. Bingham who called the 
phenomenon ‘‘pseudo-plasticity’”’ and by Wo. 
Ostwald, who called it “structural viscosity.” 

A rational solution of the problem has there- 
fore to start with a hydrodynamic theory of 
dispersed systems. 

The methods of classical hydrodynamics 
allow one, of course, to calculate in principle the 
dissipation of energy of a flowing liquid in which 
particles of given shape and size and given 
mechanical properties are suspended. Actually 
the mathematical difficulties are however very 
great and sometimes insurmountable. The vis- 
cosity of very dilute systems has been determined 
by this method in the cases where the suspended 
particles are (a) rigid spheres, (b) ellipsoids, (c) 
liquid drops. 

The method was first introduced by Einstein 
who calculated the case (a). 

It is based on the following considerations :—If 
we have a flowing liquid and we imagine part of 
the liquid to be bounded by a spherical surface, 
this sphere will be deformed through the flow 
into an ellipsoid. If the flow is laminar the main 
axis of the ellipsoid will be inclined to the 
direction of flow by an angle of 45 degrees and 
besides the deformation there will also be a 
rotation of the sphere. If however this sphere is 
not part of the liquid, but is a rigid body, it will 
be able to rotate with the liquid but will not be 
deformed. The presence of the rigid sphere 
therefore modifies the deformation of the liquid, 
i.e., it interferes with its ‘‘free’’ flow. By cal- 
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culating the dissipation of energy of a volume of 
liquid containing a number of such rigid spheres 
and comparing it with the dissipation of energy 
at ‘‘free’’ flow Einstein found that the viscosity 
Neoln, of the solution is 


Nsoln. — Nsolv. | 1 +C, 2) (1 —2C,), (58) 


where neoiy. is the viscosity of the solvent and 
C, is volume concentration of the solution. This 
equation is only valid where the particles are so 
wide apart that they have no influence one upon 
the other. As furthermore C, is very small, we 
can write 


Nsoln. = Nsolv.(1+2.5C,+5C,?+: ++) (59) 
or approximately 
Neon. = Noolv.(1+2.5Cy). (60) 


The remarkable feature of Einstein’s equation 
is this, that the size of the suspended particles is 
irrelevant as long as all are spheres, i.e., all are 
geometrically similar. 

The same method was used by G. B. Jeffery 
for ellipsoidal particles. Here, however, the 
mechanism is complicated through the rotation 
of the particles. The magnitude of interference 
of the ellipsoidal particle with the free flow 
depends, of course, on the position of the particle 
in the liquid, i.e., on the angle between the main 
axis of the ellipsoid and the direction of flow. 
This angle, however, is constantly changed 
through the rotation of the particle in the 
liquid. If the particle is a sphere all axes of the 
particle are geometrically equivalent and the 
rotation is of no influence. Jeffery did not bring 
the problem to a solution, but he found that in 
this case the viscosity can also be expressed in 
the form 


Nsoln. — Nsolv.(1+aC,), (61) 


where a is a constant which depends on the ratio 
f of the large to the small axis of the ellipsoid. 
It is of course assumed that all particles are 
geometrically similar, i.e., in the case of a rota- 
tional ellipsoid that f is a constant for every 
solution. 
Eq. (61) can also be written in the following 
form 
Nsoln. — Nsolv. Nsoln. 


Napec. _ — 


Nsolv. 





—1=n,,—1=aC,,(62) 
Nsolv. 
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where spec. is called the “‘specific’’ viscosity and 
Nrel. iS the “‘relative”’ viscosity. 

R. Eisenschitz calculated a for very long ellip- 
soids, i.e., when f is very large. He found 


ao= (1.15/m)(f/In 2f). (63) 


The calculation was carried out under the 
assumption that all values of a are equally 
probable at the beginning of the movement and 
that the particles become oriented due to the 
flow of the liquid. If, however, there is a 
Brownian movement of the particles. which 
completely preserves the isotropic distribution 
of the main axes of the particles, he found in a 
later paper 

@,,= f?/15(In 2f—3/2). 


We have denoted the coefficient a of Eq. (63) 
by ad» because it corresponds to a zero influence 
of temperature, while a@,, corresponds to an 
infinitely great influence of temperature. As can 
be seen, we have 


(64) 


a. > Qo. (65) 


Generally a will be 


ao <a <a,,, (66) 


depending on the temperature and increasing 
with increasing temperature. 


17. 


The non-Newtonian behavior of a dispersed 
system can be explained from the fact that a 
suspended particle interferes with the mobility 
of part of the dispersion medium, if it is assumed 
that this interference decreases with increased 
shear. The problem is therefore to explain why 
and in what manner an increased shear should 
decrease the interference. 

We can distinguish two sorts of interference 
which we may call (1) interference and (2) 
immobilization (after a term introduced by Wo. 
Ostwald). 

“Interference” is the mechanism which has 
been described in the preceding section. ‘‘Im- 
mobilization” concerns the fact that in ‘‘solvated”’ 
colloids part of the particle binds part of the 
dispersion medium so that a larger particle is 
formed which is suspended as a whole in the 
dispersion medium. Immobilization is therefore 
due to solvation. 
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According to Haller** we have to distinguish 
two sorts of solvation: (a) the “energetic” 
through “adsorption” and (b) the “‘steric”’ 
through geometrical inclusion. 

Lyophilic dispersions will generally show im- 
mobilization and interference as well. 

Interference can lead in four ways to a non- 
Newtonian behavior, if the particle is not a 
rigid sphere. ‘ 

(a) At rest all values of the above-mentioned 
angle a are equally probable from which there 
results a certain viscosity mo. Through flow a is 
changed, the particles become more or less 
oriented. If the orientation is complete the vis- 
cosity assumes a certain value 7,,. From the 
investigations of Jeffery and Ejisenschitz there 
follows that n,,<7o. This orientation, however, 
is counteracted by the rotational diffusion of the 
Brownian movement. Where the Brownian 
movement is negligible 7,, is established as soon 
as the steady state is reached and the system is 
therefore Newtonian with the viscosity 7... 
Where the influence of the Brownian movement 
prevails over the orientation forces yo is main- 
tained and in this case also the flow is Newtonian. 
The first of these cases is characterized by a 
vanishing stability coefficient Co, the second by 
an infinite stability coefficient. (See Fig. 8.) 
Generally there will be at every temperature and 
every shear a dynamical equilibrium between 
the forces of diffusion and orientation, with a 
corresponding 7», where o>7>»7,,. This is 
according to Staudinger the mechanism of non- 
Newtonian behavior. As has, however, been 
shown in the preceding section, the relative vis- 
cosity must in this case be dependent on the 
temperature. An increase of temperature brings 
with it increasing Brownian movement and 
therefore a shifting of 7 towards np. nrei. therefore 
increases with temperature. 

(b) The interruption of the deformation of the 
liquid produced by reaction in the particle ten- 
sional stresses. Through these stresses the 
particle is deformed. The deformation may be 
considerable. For example, it has been assumed 
in order to explain the elasticity of rubber that 
the rubber molecule is formed like a spring. In 
this case the molecule will be extended elasti- 


% W. Haller, Koll. Zeits. 57, 197 (1931). 
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Fic. 8. Fluidity as affected by the stability coefficient, C, 


cally. The tension and extension is at its highest 
when the axis of the spring has the direction of 
the principal stress. This direction is inclined by 
an angle of 45 degrees towards the direction of 
flow. As has, however, been explained before, the 
particle is not only deformed, but also rotates. 
Through rotation the spring is unloaded and 
recovers. This results in an oscillation of the 
particle, which is transmitted to the liquid and 
by this therefore the amount of dissipated energy 
is increased. If such a mechanism exists, the 
result must therefore be an increase of the 
relative viscosity, the same as through tur- 
bulence. 

(c) If the particles are threads without 
rigidity, as has been assumed by Haller the 
result is entirely different. The tension produces 
a straightening out of the thread. This decreases 
the viscosity because the particle is giving way. 
Rotation of the particle tends to restore a 
random distribution. It is difficult to judge the 
influence of change of temperature; 7,1, may be 
independent of temperature. 

(d) If the tension is very great, or the breaking 
strength of the particle very small, the particle 
breaks as soon as its axis coincides with the 
direction of main stress. The particle can also 
be broken through bending. Staudinger pictures 
his rod molecules as of very small strength and 
easily broken. Shortening of the particle decreases 


f in Eqs. (63) and (64) and therefore decreases 7. 


Whenever the fragments touch each other in a 
suitable position, the particle may be reformed. 
The probability for this is, however, small. The 
consistency curve must in this case therefore 
show a hysteresis loop. 
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TABLE II. Types of non-Newtonian behavior. 











1. Disturbance 
a. Competition between orientation and Brownian movement 
b. Elastic deformation 
c. Thread deformation 
d. Adsorption 
2. Immobilization 
a. Adsorption 
b. Steric immobilization (structural viscosity) 
a Macromolecule 
B Micelle 
y Ordered group 
6 Unordered aggregate 











(A) (B) (C) 
Variation of 7. Reversibility—absence | Existence of a lower limi 
rel. k t 
with temperature of a hysteresis loop | of non-Newtonian behavicr 

increases yes no 

? yes no 

? yes no 

nil no yes 

decreases yes yes 

? ne yes 

? no yes 

? no yes 

nil or decreases yes yes 








(2a) In the case of immobilization through 
adsorption a non-Newtonian behavior can be 
explained by a shearing off of the adsorption 
layer. Such a mechanism has been suggested by 
Hatschek and mathematically formulated by 
Reiner and Riwlin.** In this case there must be 
an influence of change of temperature. As many 
have pointed out adsorption is an exothermal 
process and should decrease with increasing tem- 
perature. In this case 71. therefore decreases, 
with increased temperature. 

(2b) The steric immobilization has been dis- 
cussed at great length by Kraemer and William- 
son. They imagine the particle to be “an 
ultramicroscopic body which is generally per- 
meated in jel-like fashion with the liquid medium. 
The union of the structural elements forming the 
particle may involve any or every type of bond 
which can operate between atoms or molecules.”’ 
The particle may be either (a) a giant molecule, 
immobilizing solvent between the branches and 
in its ring, (8) a micelle (crystallite), (7) a 
group of micelles with a regular arrangement, or 
(6) an aggregate or unordered assemblage of 
smaller particles. 

A non-Newtonian behavior can result in two 
ways: (I) through deformation, (II) through 
destruction. With regard to (I) it should be 
noted that the particle can be much deformed. It 
will therefore partially participate in the defor- 
mation of the liquid. Such a deformation of the 
particle will squeeze out of it solvent as from a 
sponge, i.e., part of the dispersion medium is 
being mobilized. With regard to (II) it should be 


noted that the deformation is accompanied by 
stresses. These stresses the particle can only 
resist to a certain limit. If this limit is exceeded 
the particle breaks and this again makes free 
part of the immobilized liquid. If the fragments 
of the particles come to touch each other in the 
course of flow, the particles may be reformed. 
“For a particular rate of shear, a dynamic 
equilibrium would be established between the 
the rate of disaggregation and the rate of reag- 
gregation.”’ It should, however, be noted, that 
the probability of reaggregation is generally not 
the same as that of disaggregation. Both are 
equal only in the case of group (2b, 8), i.e., in 
the case of unordered aggregates. Only in this 
case every meeting of the fragments in any 
position can result in a reaggregation. In the case 
of molecules, micelles and ordered aggregates a 
binding of the fragments can only take place if 
the fragments are in a certain position. Not every 
position is therefore equally favorable. If the 
consistency-measurements are made in both ways, 
i.e., the shear is increased from zero to a maxi- 
mum and then decreased to zero again, there will 
result not a single curve, but two such curves 
forming a hysteresis loop. 

It is not possible to say a priori what the in- 
fluence of a change of temperature should be in 
the case of steric immobilization. There is the 
possibility that 7,1. may not change at all with 
temperature. However, it is mostly assumed 
that in the case of an unordered aggregate the 
binding forces are so small that the aggregate will 
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partially be destroyed by an increase of tem- 
perature. 

Table II gives a summary of all these mechan- 
isms and their properties. 


18. 


The considerations of Section 17 lead to the 
conclusion that there may exist different types 
of non-Newtonian liquids. We have three criteria 
which may enable us to decide whether a given 
non-Newtonian liquid belongs to one type or 
another. They are (A) The influence of a vari- 
ation of the temperature on ‘re1.; (B) the re- 
versibility or irreversibility of a variation of 
shear; (C) the existence of a lower limit of non- 
Newtonian behavior. 

With regard to (A) it is obvious that fre1. 
cannot show a dependence on the temperature 
as long as the phenomenon of a non-Newtonian 
behavior is entirely of geometric-mechanical 
character. In this case the vibration of the vis- 
cosity of the solution with temperature is solely 
due to the variation of the viscosity of the solvent. 
(1d) represents such a mechanism, as it consists 
in a breaking, i.e., a shortening of the particle. 
Here the non-Newtonian change of viscosity is 
entirely due to a change of the geometrical shape 
of the particle. On the other hand, we can say 
that there must be a dependence on temperature 
in the case (a). A dependence must also exist 
in the case (2a). Nothing can be said a priori 
with regard to the other mechanisms. 

With regard to (B) we have already said that 
an instantaneous change of viscosity with chang- 
ing shear, i.e., without any lag is only possible 
if all positions of the particles are equally favor- 
able for a reaggregation. This is not the case if 
the particles are molecules or micelles. 

With regard to (C) it should be noted that 
wherever the non-Newtonian behavior is due to 
a partial destruction of the structure (2b), (b), 
(c), (d) or to a shortening of the primary particle 
(id, 2a) the non-Newtonian behavior must show 
a lower limit. The forces which form either the 
primary particles or the secondary structures 
can of course only be overcome when the stresses 
produced through the flow reach and exceed a 
certain value. Up to this value the behavior of 
the system is normal, i.e., Newtonian. 
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Fic. 9. Generalized Ostwald curve. 


Such a limit cannot exist in the cases (1a, b, c) 
where the smallest stress must be of influence. 

It should be furthermore noted that all 
mechanisms require that when 7,, is reached, the 
liquid behaves Newtonian. This value of n,, may, 
of course, be reached either at a finite shear or 
only at infinite shear. In both cases the con- 
sistency curve must have the shape predicted by 
Wo. Ostwald (Fig. 9). 

Experiments made up to now with non- 
Newtonian liquids have not provided enough 
information to decide which type of non-New- 
tonian behavior the liquid represents. It is 
necessary to make the observations with in- 
struments of varying dimensions in order to 
decide whether there is a wall-effect or not. It is 
furthermore necessary to use both the rotation 
type and the capillary tube type of instrument in 
order to cover a wide enough range of shear. 
Furthermore the temperature also must be 
changed. 

Such an investigation has recently been made 
with rubber-toluene solutions and it was found 
that these solutions belong to type (2b, 4). 

I want also here to thank Professor E. C. 
Bingham for the interest he has taken in this 
investigation and the valuable suggestions which 
he made in the numerous discussions. My thanks 
are also due to those benefactors and the Chem- 
ical Foundation whose financial assistance made 
possible the research of which this paper is an 
outcome. 
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Support for the assumption that rubber particles in rubber-toluene solutions are solvated 


unordered aggregates is derived from a quantitative study of viscometric measurements, which 
shows that the forces acting between the primary particles (structural elements) must be very 


small. 





FORMER communication! reported vis- 
cometric rubber-toluene 
solutions in different capillary and_ rotation- 


observations on 


viscometers. The consistency curve of solutions 
from 0.08 to 1.77 percent concentration starts at 
the origin as a straight line inclined at a definite 
angle, until it reaches a certain point a (see Fig. 
1). From there it becomes convex toward the 
stress-axis. This convex part extends to an in- 
flection point b where the curve becomes concave 
up to a point c; where it again becomes a straight 
line which by extrapolation passes through the 
origin. The solutions therefore behave as New- 
tonian liquids between points o and a, and c 
and «©. Between point a and c they are non- 
Newtonian liquids. 

It was furthermore shown that the “‘relative’’ 
consistency curve, i.e., a curve in which the 
shearing stress is reduced or the rate of shear 
increased in proportion to the viscosity of the 
solvent, is independent of the temperature at 
which the viscosity measurement is made. From 
this latter property the conclusion was drawn 
that the phenomenon of non-Newtonian behavior 
is a purely mechanical one. It cannot be due, as 
has been suggested by Staudinger,” to a com- 
petition between orientation due to laminar flow 
and rotational diffusion due to Brownian move- 
ment, because the magnitude of the latter 
depends on the temperature. 

The independence of the temperature also 
rules out the assumption of an adsorption layer 
of variable thickness, proposed by others, be- 
cause the thickness of such an adsorption-layer 
can hardly be independent of the temperature. 
It was accordingly concluded that the non- 


1M. Reiner and R. Schoenfeld-Reiner, Koll. Zeits. 65, 44 
(1933). 

2H. Staudinger, Die hochmolekulare org. Verbindungen, 
p. 189, Berlin, 1932. 


Newtonian behavior must be due to the rubber- 
particles in solution being, “‘gel-fragments,”’ or 
secondary particles ‘permeated in gel-like fashion 
with the liquid medium” as suggested by 
Kraemer and Williamson.’ This model explains 
the existence of point a. It is defined by the 
strength of the bond between the primary par- 
ticles or structural elements. Until the stresses 
produced in the secondary particles by the flow- 
deformation reach this breaking strength, the 
solution must naturally behave as a Newtonian 
liquid. When the shear is increased beyond this 
point, the secondary particles are broken up to 
form primary particles, releasing at the same 
time liquid medium. This increases the fluidity 
of the solution and causes the curve to become 
convex. The fragments formed through the 
process just described, reaggregate when coming 
in contact in the course of flow and a dynamic 
equilibrium may establish itself between the dis- 
aggregation and reaggregation. 

According to the bond, which operates between 
the primary particles, the secondary particles can 
be either macromolecules, micelles, ordered or 
unordered groups. The first three types however 
have the property in common that the fragments 
will not reunite at every chance encounter, but 
only when they meet in certain orientations to 
each other. As a result, a time-lag should be 
observed if a closed consistency curve were 
described, i.e., the consistency-curve should 
show a hysteresis loop. This was not found to be 
the case by our observations. Regardless of 
whether the rate of shear was being increased or 
decreased, the observational points fall upon a 
single smooth curve. This excludes the possibility 
of ordered aggregates and therefore leaves as the 


3E. O. Kraemer and R. V. Williamson, J. Rheology 1, 76 
(1929), 
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Fic. 1. General shape of the consistency curve. 


only possibility of unordered aggregates of 
primary particles, held together by cohesive 
forces, the primary particles themselves being 
either macromolecules, micelles, or ordered 
groups. These preliminary conclusions drawn 
from the observations were purely qualitative. 
In the present paper we will treat some of our 
experimental data quantitatively. 

The quantitative relation which is most com- 
monly looked for in viscometric observations 
with solutions of different concentrations is one 
between viscosity and concentration. When, 
however, a non-Newtonian liquid flows through 
a tube (or any other viscometer), there prevails 
in each telescopic layer a different rate of shear 
and therefore a different viscosity. It is therefore 
generally impossible to speak of the viscosity of 
the liquid at a certain efflux, corresponding to a 
certain point on the consistency curve. While 
rubber toluene solutions are generally non- 
Newtonian liquids, it was found that there exist 
two regions in which the solutions behave as 
Newtonian liquids. In these regions definite 
constant viscosities therefore exist, which we 
may denote by mp and 7n,,. 

The viscosity, 70, at rest and up to the point a 
may be called the ‘‘zero-viscosity.”’ The relative 
zero-viscosities of different solutions were deter- 
mined from our observations and are shown in 
Table I. Table I also shows the “‘specific’’ vis- 
cosity, which is defined by 


specific »=relative 7—1 (1) 
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TABLE I. Relative and specific zero-viscosities and fluidities. 








No. of 
solution II 19 


~~ 


6 4 3 2.1 2 





Concentration 


10? cy 1.71 1.48 0.98 0.97 0.23 0.17 O.11 0.08 
rel Ny 185 108 36.2 27.3 2.55 2.15 1.76 1.60 
rel %o 00541 .00943 .0276 .0366 .3922 .4651 .5681 .625 
Spec ny 184 107 35.2 263 1.55 1.15 0.76 0.60 
a, 10800 7240 3600 2700 670 680 690 750 











and a, the specific viscosity per unit of concen- 
tration, i.e., 


a=specific 7/C», 


(2) 


where c, is the volume-concentration of the 
solution. 

In Fig. 2 rel no and a» have been plotted as 
functions of c,. As can be seen from this figure, 
rel no is linear and ap constant up to c,=0.3 
percent, corresponding to a value of spec m)=2. 
This is remarkably high, as it is usually assumed 
that with solutions of high molecular substances 
the simple linear relationship has its limit with 
spec no=0.3.4 The constant value of ao is 700. 
This can be interpreted as follows: (a) For very 
dilute suspensions of spherical particles Einstein 
derived the formula 


specific 7 =2.5C,, 


(3) 


where C, is the concentration of the suspended 
particles. If the suspended particles are solvated, 
i.e., if they immobilize part of the dispersion- 
medium, we have 


C,=be,, (4) 


where b>1. b has been called the specific hydro- 
dynamic volume of the material. 

Introducing the value for C, from (4) into (3) 
we get 


specific 7 =2.5bc, (5) 
and by considering (2) 
a=2.5b or b=a/2.5. (6) 
In our case we would have 
bo = 700/2.5 = 280, (7) 


i.e., if the dispersed solvated rubber particles 
were of a spherical shape, the volume of the 
immobilized liquid would be 279 times the 


4E. Guth and H. Mark, Ergebnisse, p. 118, Berlin, 1933. 
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Fic. 2. Relative zero viscosities for different concentrations. 


volume of the dispersed rubber and the particle 
concentration of our solutions, 2, 2.1, 3, and 4 
would be 22, 31, 48, and 65 percent, respectively. 
This however is not possible. Einstein’s equation 
is valid for very dilute solutions only, and it is 
well known that the viscosity increases linearly 
with concentration at small concentrations only. 
As has been said before the increase of viscosity 
is linear and quite small up to c,=0.30 percent. 

If the real concentration corresponding to 
c, =0.30 percent was 2800.3 =84 percent, this 
would be out of question. The actual ratio of 
volume of immobilized liquid to volume of 
rubber must therefore be much smaller than 279. 
(b) Einstein’s® formula (3) can be generalized to 


specific 7 =AC, (8) 


for particles of general (non-spherical) shape, 
where A is a different constant for every dif- 
ferent shape of the particle, being equal to 2.5 
in the case of a sphere. Eisenschitz® could show 
that under certain assumptions 


1.15 / 1 
A =-—.---——— (9) 
x d \n2l/d 


5A. Einstein, Ann. d, Physik 19, 289 (1906); 34, 591 
(1911). 
® R. Eisenschitz, Zeits. f. physik. Chemie A158, 78 (1931). 
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in the case of very long rotational ellipsoids of 
length / and thickness d. 

If c, were the true volume-concentration of 
the dispersed phase, i.e., if the dispersed rubber 
was not solvated, and if we had reasons to 
assume that the particles are long rotational 
ellipsoids, we could calculate from Eq. (9) that 





1 1 7007 
“ve = = 1910 (10) 
d\nI/d 1.15 
or 
1/d = 20,000, (11) 


which means that the particles were of a length 
of 20,000 times their thickness. We know, 
however, that this cannot be the case, as the 
particles are solvated. We have therefore to 
picture the particles as solvated to less than 279 
times and elongated to less than 20,000 di- 
ameters. 

Table II shows possible combinations of 
solvation and elongation. Under C, the particle 


TABLE II. Possible combinations of solvation and elongation 








for a=700. 
l 280 100 10 1 
l/d 1 20 1500 20,000 
‘ai 84% 30% 3% 0.3% 








concentration up to which the viscosity-concen- 
tration relation is linear, is shown. Beyond this, 
we cannot at this stage, draw any conclusions 
about the shape of the particle and the quantity 


of liquid imbibed. 


This is how far the information provided by 
the first straight part of the consistency-curve 
up to point a carries us. Next the establishment 
of a relation between the position of the point a 
and the concentration may be undertaken. 

We said that the existence of a point a (Fig. 1) 
must be interpreted as indicating the starting of 
a mechanical disaggregation of the suspended 
rubber-particles, through the breaking off of 
primary particles. This breaking off is caused by 
the stresses, produced in the particle through 
the flow of the solvent, when these stresses 
exceed the breaking strength 8 of the particles. 
We know from our data (see Fig. 3) the shearing 
stress (P) and relative rate of shear (Vye1) at 
which a is reached in the case of some of the 
solutions and it should be possible to find a 
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Fic. 3. Consistency curves at low rates of shear. 


relation between these magnitudes and the 
breaking strength of the particle. Such a relation 
can be theoretically derived on the basis of 
Einstein’s calculation of the viscosity of a solu- 
tion as shown in the Appendix. The result is as 
follows: Let V;: be the relative shearing velocity 
of laminar flow, then a tensional stress S; is 
developed in the particle, the magnitude of 
which is 

S:=1V rei, 


(12) 


where 7 is a numerical factor, which in the case 
of spherical particles is equal to 5 and for 
“elongated” particles is greater than 5. The 
stress acts in a direction inclined under 45° 
against the direction of flow. It is independent of 
the concentration, the size of the particle and is 
the same in every point of the particle. 

The particle begins to disaggregate as soon as 
S, reaches the breaking strength 6. Disaggre- 
gation will therefore start at the limiting relative 
shearing-velocity 


Vret= p/n. 


The limiting relative shearing-velocity is there- 
fore independent of the concentration. 

On the other hand P, the limiting shearing 
stress, is 


(13) 


P= Nrel* V ret= (p/n) nret (14) 


and as rei depends on c,, P depends in the same 
way on c,. This seems not to be confirmed by 
Fig. 3. Here it looks as if on the contrary P was 
independent of c, and V,., was decreasing ac- 
cordingly. It should, however, be born in mind 
that the consistency-curves (Fig. 3) refer to the 
solutions II and 19 with the high concentration 
of 1.71 and 1.48 percent. For such concentrations 
our Eqs. (13) and (14) are not valid. Unfor- 
tunately, in the more dilute solutions the point 
a could not be reached at the maximum shear 
produced in the Couette-apparatus. 

Solutions II/2 and II/18, both with the same 
concentration of 0.08 percent were however 
measured in the Couette and Bingham-viscom- 
eter and it was therefore thought possible to 
find P and V,.; by a graphical comparison of the 
measurements in both instruments. Such a 
graphical comparison is made difficult by the 
fact that the shear produced in the capillary- 
instrument is so much larger than the shear 
produced in the rotation-instrument that a 
graph, as usually plotted becomes unwieldy. 
This difficulty is partly overcome by plotting on 
a logarithmic scale. In the same way the meas- 
urements of solution III/14 and 13 and II/15 
in the Bingham-instrument were compared with 
the hypothetical between the measurements in a 
Couette-instrument derived by intrapolation 
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Fic. 4, Limiting rates of shear and shearing stresses. 


between the measurements made with solutions 
I1/17, 11/6, I1/4 and 3. The results are shown in 
Table III and in Fig. 4. Fig. 4 confirms our Eqs. 


TABLE III.* 











No. of 
solution II 19 14 13 15 18 
102 cy 1.71 1.48 0.39 0.38 0.21 0.08 
Vrel 0.034 0.058 2.55 7.79 7.35 
P 6.30 6.30 13.5 22.7 11.0 








* Because of the method used in the determination of Vrel and P, the 
values given for solutions 14, 13, 15 and 18 are very approximate. 


(13) and (14) for very dilute solutions. At con- 
centrations over 0.25 percent, however, both Vre1 
and P fall rapidly, the latter tending towards a 
constant. The concentration of 0.25 percent is 
also approximately the concentration at which 
Vet ceases to increase linearly (see Fig. 2), i.e., 
this is the concentration at which the suspended 
particles influence one another in their hydro- 
dynamic behavior. It may be expected that such 
an influence should bring with it a rapid decrease 
of V,.; and P. 

The results of the proceeding section confirm 
our interpretation of the point a of the consist- 
ency curve and therefore also of our model for 
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the suspended particle and the mechanism of 
the non-Newtonian behavior. For very dilute 
solutions V;.; tends to become constant with 
the value 7.70. According to equation (13) we 
have 


8 =7.70n dynes /cm?. 


We do, however, not know 2, as we do not know 
the shape of the particle and are therefore not in 
a position to calculate 8, the breaking strength 
of the suspended particles. This breaking 
strength is not identical with the cohesion of the 
rubber because in the solution the cohesion 
forces in the rubber are counteracted by the 
adhesion between rubber and solvent. Moreover 
the inclusion of the solvent reduces the massive 
section of the particle. It is not possible to cal- 
culate this reduction, as we do not know the 
volume of the solvent, which is immobilized by 
the rubber particles. We can nevertheless say 
that the bond between the primary particles, as 
measured by 8, is very small. If we take n=5 we 
get 8=38 dynes/cm?’ and this is about 10-7 the 
tensional strength of rubber. It can hardly be 
assumed that » should be so large as to make up 
this difference. We must therefore conclude that 
the forces which combine the primary particles 
to the secondary particles are very small. 

The present quantitative analysis of the 
results of viscometric measurements with rubber- 
toluene solution has furnished new support for 


as 


Fic. 5. Stresses in a spherical particle. 
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the model of solvated unordered aggregates, sug- some other measurements with the viscometric 
gested for the suspended rubber particles. It was measurements as, e.g., diffusion measurements. 
however not possible to calculate the amount of I am glad to express my gratitude to Dr. E. C. 
immobilized solvent and therefore also not the Bingham for the interest taken in this research 
size of the suspended particles. It is thought that and to those benefactors who have made it 
this can only be accomplished by combining _ possible. 





APPENDIX 


To calculate the stresses produced in the particle through its interference with the free flow of the 
liquid, we start from the following equations of Einstein’s paper (page 295) 


’ * 
: Ou dv dw X,¥,2 
X:= p—2k—; Y;=Z,= —_ (—+—) é, n; i? ’ (1) 
dg ag On ) 
u, VU, W 
where u, v and w are according to Eq. E(6)** 
Fe 5 P* | A, B,C 


u= At—~ —¢(AE+By?+Ct*)+-—(AP+Bytt+ce)—-—4t (tans) (2) 
2 p° 2 p? p® 

















and 
p= (B+y%+ 59) (3) 
From (2) there follows 
Ou 5 P8 5 P3 5 Op 5 P5 
=A-—- —s 24§-- — (Ae +Br?+Ce +s PreAe + Bry? + Cg") — ie —§-2Aé 
ae 2 p® 3 p° p° dE 2 9’ 
§ P* 7 0p P*® 5 dp 
a Ae +Brt+ Ci) —*P> g(Ae? + Br? +Cy?)—— —-—A+P*At——. (4) 
p? p> de p® p° a 
From (3) there follows 
dp/dE=2E/2(E+7°+5)*= E/p. (S) 
Introducing 0p/0é from (5) and putting p=P, Eq. (4) becomes 
Ou 5&? Lfiensibes 
‘at -—"|4-<ar +Bn? +c7)| g, UE c . (6) 
Slop A, B,C 
Similarly we get 
Ou 5nd 1 veins 
~ =] 4-4 +B +c] & mo} (7) 
aad A, B,C 
Furthermore we have 
1/p=1/(&+7?+$7)! (8) 
and therefore 
d(1/p) E O*(1/p) 2&?—f5? 3£%—p? 
=-—; — = . (&, , f). (9) 
d& p® 0g p° p° 


If we introduce the expressions for 0?(1/p)/d#, 0°(1/p)/dn’, and 6?(1/p)/d¢? into the first of Eqs. 
E(5) we get 


i. ay - be permuted cyclically in order to obtain two more 
*1£, », € } means that X, Y, Z, &, n, ¢, u,v, wshould analogous equations. 
u, v, Ww 


** FE (6) means Eq. (6) in Einstein’s paper. 
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5 f 3¢?— p? 3n?— p? 3¢?— p? 
p= —-kP3\ A +B +C (10) 
3 | p® p® p® 





and considering Eq. E(2) and putting p=P this gives 


Sk 
= — + By’? +Cf*). (11) 





p=P 


We now introduce the values from (6), (7) and (11) into Eqs. (1) and get the following expressions for 
the components of the stresses acting on the surface of a small spherical particle 





5k e 
X;= -|(: -2=) (AP + Bat +-CH*) +248] Xx, Y,Z2 
(: ”, 2) (12) 
A, B,C 


Snckp 2 
¥;@ 2," + —(AR+Br+ Cr) +4 
P2 | P? 


In Eqs. (12) we have gained expressions for the six components of the stress tensor in every point 
of the surface of the particle with the coordinates £, n, ¢, where 


SP-++-y?-++- fem P*, (13) 


We now want to know the three components of the forces per unit area of the surface of the particle, 
parallel to the axes £, 7 and ¢ and exerted by the liquid upon the particle. 
They can be calculated from the Eqs. (12) by means of Einstein’s equation page 295 


xX, Y,Z 
). (14) 


X,= —(X:é/P+X yn/P+X¢f/P) ( 
We get 


5k 


é 
~X.=— (AP+ Bat + Cx E22 +r +# |+e(-24e+30+-00)| 


Ské . 


: SRE 
= pe art £7)+C(n?+ Jarl AP (+ tA +B+C)] (15) 


and if we consider Eq. E(2) this gives 


X n= 5kAE/P, Y,=5kB,/P, Zun=5kC¢/P. (16) 
For every value of — between — and =P we have a calote on which the hydrodynamic force 
P 
X= ( 2(P?—£)'tX ads (17) 


acts (see Fig. 5). It can easily be seen that 
ds=([P/(P?—é)! ]dé. (18) 
Introducing ds from (18) into (17) we get 


, SRA 
X= 2Prf X ,d§=2Pr 
g 





f tdt= 5xkA(P?— £2), (19) 
g 
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This force produces internal stresses in the particle which in a section normal to é axis and at the 
distance ~ are equal to 


& 2 ¢ 
=X/(P?—#)r=5kA ( ). (20) 
Pi ( A.B.C 


This result is quite remarkable. Eq. (20) says that the stress is independent of P and &, i.e., the 
reaction of the liquid against the interference of its free flow by a spherical particle, produces in this 
particle internal stresses which are independent of the size of the particle and are the same for every 
parallel section through the particle. 

We have derived the stresses produced in the spherical particle under the assumption that there 
is only one particle suspended in the liquid. We now have to proceed to the case when there are 
many particles suspended in the liquid, which case is dealt with in §2 of Einstein’s paper. In equations 
(20) A, B, C are the main extension-velocities of the liquid, if there was no spherical particle sus- 
pended. These have to be replaced by the main extension-velocities A*, B* and C*, of the “‘solution.”’ 
We get accordingly 


‘ & mn ¢ 
pet=5kA* ( ). (21) 
A*, B*. C* 
where the symbol* refers to the solution. According to our definition of V1, we have 
A rei=kA* (A, B, C), (22) 
and therefore 
pt=5A (; n, ') (23) 
=5Are . 
A, B, C 


From considerations of mechanical similarity we may conclude that for particles of any shape Eq. 
(23) can be generalized to 
S:=nV ret; (24) 


where S;, is the tensional stress produced in the particle, V;.1 is the relative shearing-velocity and n 
is a numerical factor which in the case of a sphere is equal to 5 and for “‘elongated”’ shapes must be 
greater than 5. 
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INTRODUCTION 


HEN making rheological measurements on 

non-Newtonian liquids, especially those 
that are thixotropic, it is highly desirable to use 
an instrument in which the liquid is subjected to 
continuous and uniform or essentially uniform 
stresses. Otherwise the data cannot be inter- 
preted without ambiguity or questionable as- 
sumptions. The rotating cylinder viscometer in 
its various forms is the only viscometer so far 
developed that meets this requirement; but it is 
usually made difficult to clean and is complicated 
mechanically by guard rings placed at the ends 
of the inner cylinder to eliminate end effects. The 
upper guard ring can easily be dispensed with if 
the liquid level in the viscometer is kept below 
the top of the cylinders. The lower guard ring, 
also, has recently been eliminated by an in- 
genious design, used by Mercier,! in which an 
air pocket in a concavity in the lower end of the 
inner cylinder practically eliminates shearing 
stresses over the bottom surface. 

The essential feature of a new design described 
in the present article consists in such a formation 
of the lower ends of the two cylinders that the 
ratio of the clearance between them at any point 
to the radial distance of the point from the axis 
is constant and is equal, at least approximately, 
to the ratio of the clearance to the radial distance 
at the sides of the cylinders. The rate of shear in 
the liquid at the bottom will therefore be essen- 
tially the same as at the side. The outstanding 
feature of the rotating cylinder viscometer is 
thereby retained, while the guard rings are 
eliminated without the introduction of an air 
pocket at the bottom. 

In the form of the instrument described below, 
called the conicylindrical viscometer, the bottom 
surfaces of the two cylinders are cones of slightly 
different angle. The hydrodynamical analysis of 
the conical region will first be developed, after 
which the instrument, as built, will be described 
in detail. 


'P, Mercier, J. Rheology 3, 391 (1932). 


HYDRODYNAMIC ANALYSIS 


The complete notation is given in the table at 
the end of the paper. Referring to Fig. 1, let a 
and b be the radii of the inner and outer cylinders, 
respectively, and a and §6 the corresponding 
angles of the conical bottoms. We need to analyze 
the hydrodynamics of the region between the 
conical surfaces only, the solution for the 
cylindrical surfaces being known. Let 7, 0, ® be 
a system of spherical coordinates as indicated. 

Suppose that the lower cone is rotating in the 
positive direction with angular velocity w, that 
the upper cone is stationary, and that the space 
between them is filled with a liquid of fluidity ¢. 
The fluidity is not assumed to be independent of 
shearing stress. 

It can be shown from the complete system of 
stress equations and boundary conditions that, 
inertia being neglected, the liquid at any par- 
ticular point is rotating about the axis of the 
system with a velocity, v, which is a function of 
the coordinates r and @; that the other two 
velocity components are zero; and that the only 
shearing stress different from zero is the one 
acting in the #-direction on a 6-surface. All these 
results are as would be anticipated; and they 
will here be assumed without proof. 

The equations? which then remain to be 
solved are 


dp/d@+2p cot 6=0, 


| (1) 
(1/r)(dv/d0—v cot 0)= pg, 


where p=the shearing stress in the #-direction 
on a 6-surface and v=the circumferential, or 
’-velocity of the liquid. 

For the purpose of the analysis it is assumed 
that the cones are of infinite radius. The re- 
maining boundary conditions are: 


2 See, for example, Love, Mathematical Theory of Elas- 
ticity, p. 89, 3rd edition, Cambridge University Press. Eq. 
(1.1) is the form taken by the equation of equilibrium when 
there is no acceleration and other conditions are as have 
been assumed here. The left member of Eq. (1.2) is the rate 
of shear in the medium. So also is the right member, by 
definition of the fluidity. Cf. Love, p. 56. Love’s equations 
remain unaltered if we interpret strains as rates of strain. 
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v(r, a)=0, 
: (2) 


v(r, B)=rw sin B. 


The solution of Eqs. (1) subject to the bound- 
ary conditions (2) is 














{ (vd6/sin* 6) 
v=rw sin 0 —] , 
f (gd6/sin* 6) 3) 
= B . ‘ 
sin? af (od@/sin*® @) 


We now wish to determine a and 8, or their 
difference, so as to satisfy the condition that the 
shearing stress and rate of shear between the 
conical surfaces shall be the same or approxi- 
mately the same as between the cylindrical 
surfaces of the viscometer. To be more specific, 
let us seek to satisfy the condition that the 
shearing stress between the cones at the mean 
value of @ shall equal S, the mean shearing stress 
between the cylinders, or 


po= 5S. (4) 


In an analysis of the cylindrical viscometer* 
it is shown that, except for higher terms in « 
which may be neglected in the present analysis, 


*M. Mooney, J. Rheology 2, 210 (1931). 
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the fluidity at the mean stress, go, is given by 
yo w eS, (S) 


in which e=(b?—a’)/(b?+<a’). 
Then from Eqs. (3.2), (4) and (5) we have 


® 








w 
p = Edo. (6) 
sin? a | (ed0/sin® @) 


Furthermore, as a consequence of (4) we must 
have (09) =¢o. If € is small, as it is in a well- 
designed viscometer, the angle (8—a) will like- 
wise be small. Hence, by Eq. (3.2), the range in 
p will be small. Therefore, even if ¢ is not con- 
stant, the range in ¢ will be small unless there 
is a very unusual discontinuity in the fluidity- 
stress curve of the liquid. As a valid approxi- 
mation, the ¢ under the integral in Eq. (6) may 
therefore be set equal to the constant go. We 
then obtain: 





. 8 dé 
e= sin? oof . (7) 
a sin’ 6 
If we now set: 
6/2= B—0)= %— a, (8) 


the integral in Eq. (7) can be developed as a 
power series in 6 (expressed in radians), and by 
reversion we finally obtain: 


5=e sin 69+ €3(2 sin? 09—3 sin 0) +--+. (9) 


6) being chosen arbitrarily, Eq. (9) then 
determines 6, the angular separation of the 
cones, so as to satisfy the condition that the 
shearing stress and rate of shear should have, to 
the first approximation, the same values in the 
conical and in the cylindrical regions of the 
viscometer. 

Our next step is to calculate the torque re- 
sulting from the integrated shearing stress acting 
on a conical surface in the viscometer. If we 
consider the situation at the corner, where the 
cones and cylinders meet, we conclude that the 
actual torque would be slightly greater than 
what we would calculate for the upper cone and 
slightly less than what we would calculate for 
the lower cone. A similar statement applies to 
the sides of the two cylinders, which are of 
slightly different length. Therefore, without 
attempting to analyze further the corner effects, 
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we calculate the torque acting across the surfaces 
where the stress has the mean value S. In the 
cylindrical region, this surface has the radius 79, 
given by: 

ro? = 2a7b?/(a?+-5?). (10) 
In the conical region the surface is defined by 
6=6,)=(a+)/2, as previously indicated. 

The torque on the bottom is then given by: 


ro/sin 00 
T,= f 2rr*po sin? Oodr. (11) 
0 


Integrating and using Eqs. (4) and (5), we obtain, 


4rwr,a*b? 
T= - ° 
3g sin 09(b?—a?) 





(12) 


The torque acting on the cylindrical surface is, 
in our present notation, 


4rlowa?b? 
T,.=———__, (13) 
¢o(b? —a*) 


where /, is the effective length of the cylindrical 
region, measured from the upper liquid level to 
the intersection of the cylinder of radius ro with 
the cone 6=46. Jo, the immersed length of the 
inner cylinder, is given in terms of ],; by the 
relationship 


lo—l,=a cot a—?ro cot I, (14) 


which can be reduced by approximations to the 
form 
b—a 
b= 1,+——_(1 — COs 60). 
2 sin 4 


(15) 

By adding the members of Eqs. (12) and (13) 
an equation for the total torque, 7, is obtained, 
which can be solved for go, the mean fluidity, 


giving: 
w Anlya*b? ro 
eon (14 ). 
T b*?-—a? 3lo sin 4 


An explicit formula for the mean shearing 
stress, S, can likewise be obtained by suitable 
algebraic manipulation of the previous equations. 


(16) 


T 


S= : 
2rrolo(1 +7o/3lo sin 60) 





(17) 
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DESCRIPTION OF THE INSTRUMENT 


A conicylindrical viscometer is shown in Fig. 2, 
built primarily for use with rubber latex. The 
outer cylinder, 1, is the stator. The inner cylinder, 
2, is the rotor, and is driven by a pair of equal 
weights producing a torque in the same direction. 
Only one of the weights, 3, is shown in the 
drawing, suspended by a silk thread over the 
pulley 4. The stator is surrounded by an asbestos- 
lagged water jacket, 5. The jacket is provided 
with a thermometer hole and inlet and outlet 
for circulating water, which are not shown. The 
rotor turns on a jewel bearing, 6, at the bottom 
and has a steel conical bearing, 7, at the top. 
The trough, 8, surrounding the stator at the top, 
is to catch the overflow liquid from the viscom- 
eter when the rotor is lowered into place. The 
cover, 9, formed in two halves, serves to prevent 
evaporation and skin formation on the upper 
surface of the liquid. 

The radii of the cylindrical sections of rotor 
and stator are, respectively, 2.1 cm and 2.0 cm, 
leaving a clearance of 1 mm. The bottom of the 
stator is flat. The bottom of the rotor is a cone of 
half-angle 87° 2’, as required by Eq. (9). The 
inside corner of the stator where bottom and 
side meet is rounded slightly to facilitate clean- 
ing. Since the clearance between the bottom 
surface of the rotor and stator is important, a 
metal ring is provided which has the form of the 
space between the two bottom surfaces; and 
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Fic. 2. Details of the conicylindrical viscometer. 
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the height of the lower bearing is adjusted by 
means of a screw while the ring spacer is in place. 
Metal parts of the viscometer are made of 
stainless steel. 

All moving parts of the viscometer are made 
as light as possible, consistent with appropriate 
strength. The rotor is hollow and weighs 45 g 
when immersed in water up to the working level. 
The pulleys are of fiber and turn on steel conical 
bearings. The most satisfactory thread found for 
supporting the driving weights is No. 60 silk 
tatting thread. 


CALIBRATION 


Three instruments of the design shown in 
Fig. 2 were calibrated at the same time with the 
same series of glycerine solutions. These solutions 
were made up from Armour and Company 
“Chemically Pure,” U.S.P. glycerine and the 
specific gravities determined at 25° by pyc- 
nometer. The percentage glycerine and viscosities 
were then obtained by interpolation from 
Sheeley’s tables.‘ 

The temperature in the water jacket of the 
viscometer during calibration was determined 
by a thermometer graduated to 0.2°C, which 
was checked against a Bureau of Standards 
thermometer. The instrument constant k was 
defined by the equation: 


n=1/go=kT/w, (18) 


in which 7 is the mean viscosity. 

Eliminating go7/w between Eqs. (16) and 
(18), we express & theoretically in terms of the 
dimensions of the instrument: 


b?—a? 


k= ’ 
Anlya?b?(1 +7 o/3lo sin 60) 





(19) 


Table I, in columns 5, 6 and 7, gives the experi- 
mental values for k calculated by Eq. (18) from 
data obtained with a series of different glycerine 
solutions, and also, in the last row, the theoretical 
values calculated by Eq. (19). The theoretical 
values are slightly different for different instru- 
ments because of small differences in the dimen- 
sions of the finished instruments. Each experi- 
mental value of k shown is the average of either 


4M. L. Sheeley, Ind. Eng. Chem. 24, 1060 (1932). 
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TABLE I. Experimental and theoretical values of instrument 
constant k. 








: Values of kX104 
Viscosit y-Poises Instrument 


Solution Glycerine 20.0°C 22.5°C No. 1No.2 No. 3 





A 94.65 5.03 4.18 1.79 1.98 1.81 
B 92.37 3.50 2.88 1.86 1.95 1.86 
§ 90.24 2.44 2.02 1.96 2.06 2.02 
D 80.50 0.656 0.557 1.92 1.98 1.94 
E 69.12 0.213 0.189 1.89 1.90 1.88 
Experimental average omitting A 1.91 1.96 1,92 
Theoretical 1.93 1.94 1,90 








three or four values obtained with different 
driving torques. The variation of these individual 
values from the mean was about 0.5 percent. 
The fact that the checks obtained on the same 
solution with different instruments are better 
than those obtained with different solutions on 
the same instrument indicates that the vari- 
ations are largely due to errors in the assumed 
viscosities of the solutions rather than in the 
data obtained with the viscometers. 

In the investigation of the friction in the 
instrument when in actual operation, the most 
reliable test was given by the stress-rate of shear 
curve, shown in Fig. 3, when the instrument was 
filled with water. The stress intercept shown is 
equivalent to a driving load of about 0.2 gram. 
It will be noticed that at higher loads turbulent 
flow sets in. For this reason it was not possible 
to obtain reliable determinations of the friction 
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Fic. 3, Measurements with water in the conicylindrical 
viscometer. 
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With the two driving loads 
opposing each other and with only a drop of 
liquid in the viscometer on the lower bearing, 
the friction was 2.3 grams at zero load and in- 
creased about 0.5 g per 100 g total load. The 
actual friction in the instrument while in use 


at high loads. 


would be less than this because of the bouyant 
action of the liquid on the rotor and because the 
driving thread tends to bind more when the 
loads are opposing each other. 

The conicylindrical viscometer in its present 
form is thus shown to be an absolute instrument 
of considerable precision in measuring viscosities 
of 0.05 poise or more; and in this range it has 
also been found more convenient to use than the 
capillary viscometer previously used for the 
same materials. 


SUMMARY 


A new rotating cylinder viscometer is described 
in which the bottom surfaces of the two cylinders 
are cones of slightly different angle. The hydro- 
dynamical theory of the conical region is de- 
veloped, and a formula is given for the proper 
angular separation of the conical surfaces to 
produce the same rate of shear between the 
cones as between the cylindrical sides. The 
theoretical relation between viscosity, applied 
torque and rate of rotation is checked experi- 
mentally by measurements with glycerine solu- 
tions of known viscosities. 


NOTATION 


a =radius of inner cylinder, 

b=radius of outer cylinder, 

g = gravitational acceleration, 
?=instrument constant, 
=nw/T, 

/,; =immersed length of inner cylinder, 
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lo=effective mean length of inner and outer cylinders, 
_b—a 1— cos 6 
oe sin 4 
pb =shearing stress in the conical region, 
po=shearing stress at 0=4o, 
r =radius vector, 
ro = effective mean radius of inner and outer cylinders, 
= [2a°b?/(a?+5*) }', 
S=mean shearing stress at the inner and outer cylinders, 
=T,/2zlor¢’, 
T =total torque, 
=7,+Th, 
T, =torque on side (cylindrical) surface, 
7, =torque on bottom (conical) surface, 
v=velocity of the liquid. 
a=angle of the inner cone, 
8=angle of the outer cone, 
5=68—a (in radians), 
e= (bh? —a*)/(b? +a"), 
n=mean viscosity of the liquid, 
6=colatitude angle of spherical coordinate system, 
6)=mean angle of inner and outer cones, 
=(a+8)/2, 
¢g = fluidity of the liquid, 
rate of shear 


shearing stress 
¢go= fluidity at 02=@, or at r=Po, 
&=longitude angle of spherical coordinate system, 





w =angular velocity of rotor (in radians/sec.). 


GLOSSARY OF IMPORTANT EQUATIONS 








w 4nl,a*b? Yo 

nl ‘vet ) (16) 
T b?-a? 3lo sin 4 
 d 
S= (17) 
To 

2rrolol 1-+-———— 

3ho sin Ao 
b—a 1—cos 4 
lop=1,+ . ’ (15) 





2 sin Ao 


5=e sin 09+ *(2 sin® 69—3 sin 4). (9) 
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The Flow Characteristics of Sewage Sludge and Other Thick Materials* 


WILHELM MERKEL, Stuttgart 
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OR the operation of clarifying plants, there 

was need for data on the pressure loss of 
sewage sludge in pipes, to form a basis for cal- 
culating the economical pumping of sludge to 
long distances. 

The chief difficulty lay in the greatly changing 
composition of such a sludge and the different 
effects which the separate ingredients exert on 
its flow characteristics. Therefore a purely 
mechanical investigation, no matter how care- 
fully and thoroughly undertaken, as for example 
in Chicago, Baltimore and Syracuse,” could not 
lead to a satisfactory result. The author sought 
to determine the separate relations and corre- 
lations by numerous observations, and thus came 
to results which are of significance not only for 
the sanitary engineer but also for rheologic 
research and for industry, insofar as they have 
to do with the mixing and pumping of thick 
materials. 

The tests were mostly carried out at the 
Stuttgart clarifying plant, and also partly at 
Nuremberg. The best results were obtained by 
flow tests in which sludge from the fermentation 
chamber of an Imhoff tank was discharged into 
a well until it reached the level of the surface 
of the settling basin (Fig. 1). After applying 
various corrections, useful curves were obtained 
which showed the relation between the average 
velocity of flow in the pipe and the loss of head, 
from the highest value at the start of the test 
(velocity v=1-—3 m/sec., lost head J = 6 percent) 
until equilibrium was reached (v=0, J=0). The 
special losses (at entrance, elbows, etc.) could 
be eliminated by suitable arrangement and 
determined separately. Other tests were made on 
straight pipe (20-40 m long) with constant head 
and different diameters; the distribution of 
velocities on a section could be observed in open 
channels, in straight stretches, on curves, at con- 
tractions and enlargements of section. 


* Extract from a paper by the author, under the same 
title, which was published by R. Oldenbourg, Munich, as 
Series II, No. 14, Beihefte zum Gesundheits-Ingenieur. 
Translation by W. H. Herschel. 


Finally measurements with the Kaempf vis- 
cometer®: 7 at any concentration and temperature 
made possible the final explanation of relation- 
ships which could not be determined by the large 
scale tests. 

Only the rheologically important results will 
be given here. Sewage sludge is a mixture of 
water, solutions and colloids, with coarse material 
up to a few centimeters in diameter according to 
the free width in the racks before the clarifying 
plant. Of the solid parts of fresh sludge, 20-30 
percent are of inorganic nature (fine sand and 
ashes; coarse sand is mostly separated out in a 
sand trap before the clarifying plant). The rest 
is organic refuse from housekeeping and factories, 
as well as fecal matter. The organic constituents 
have the characteristic of adsorbing a consider- 
able amount of water, i.e., micelles are formed 
with liquid envelopes. This causes the sludge to 
become viscid and slimy. During the fermen- 
tation which takes place under the action of 
bacteria, the colloids are first built up with their 
great surfaces, by which a large amount of water 
is set free. The coarse lumps fall apart into small 
pieces, go partly into solution or into colloidal 
form, and then likewise disintegrate. The organic 
material, during the fermentation, undergoes 
not only a mechanical and chemical transforma- 
tion, but also a loss of material, due to the 
release of water and of gases (mainly COs and 
CH,; all together about 5-800 liters of gas per 
kg of organic material). By the decrease in the 
organic part of the solid material of the sludge, 
the viscosity is decreased, but this increases the 
effect of the inorganic, mainly granular, con- 
stituents. As the water content of such a sludge 
decreases, it finally changes to an earthy mass. 

The derivation of an equation of flow for 
sludge from the limiting conditions is impossible, 
for such an equation cannot be satisfactorily 
derived either for water (apart from laminar 
flow) nor for dry earth. 

In the investigation of the flow characteristics 
of plastic materials, four ingredients must be 
considered, each of which obeys different laws. 
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Fic, 1. Arrangement of experimental plant. 


1. The free dispersion medium; it is in general 
a simple viscous liquid with a characteristic re- 
sistance, its viscosity yu. 

2. The bound dispersion medium, forming 
with the lyophile kernel which it covers, a slimy 
mass; here there are added adsorption and 
capillary forces which resist shearing and cause 
an increase in viscosity. 

3. Small lyophobe particles which hardly 
affect the fluidity of the system when present in 
small amounts, but only when numerous enough 
for the separate particles to touch and rub 
against one another, they cause a resistance 3 
which does not vanish even at the highest 
shearing stresses. 

4. Coarse objects up to a few centimeters in 
diameter which take up a considerable part of 
the pipe section and prevent the formation of a 
normal velocity profile. 

If we take up the quantitative consideration 
of the flow characteristics we must first be clear 
about the effect of the fourth constituent. In 
fermented sludge all the large particles of the 
fourth constituent occur only in small amounts, 
and normally in fresh sludge there is so much 
water present that here also they do not play an 
important part. We can therefore leave them out 
of consideration, or include them in the very 
large tolerance required from the variations in 
composition of the sludge. 

The action of the free water remains fixed; it 
follows Newton’s law r=udv/dr. With con- 
stituent 3 it is to be observed that coarse lyo- 
phobe particles, or those which do not adsorb 
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water, cannot be regularly distributed and 
carried along in laminar flow, since they sink on 
account of their greater weight, and at most are 
rolled along the bottom of the pipe. Only in a 
medium with structural viscosity can lyophobe 
particles be carried along at any elevation. 
However we also have the so-called ‘“‘pseudo- 
laminar” flow. As above stated, a noticeable 
effect on the flow characteristics takes place at 
a concentration such that capillary effects are 
possible in consequence of the small distances 
between particles, or the solid particles rub on 
one another. This friction is nearly independent 
of the velocity (but not of the pressure) i.e., 
7=constant. To be sure, when the solid particles 
exceed a certain size, additional resistances are 
caused from the fact that the particles are not 
spherical and become wedged in the spaces 
between other particles; this is the more so 
the longer time they have, or the less the 
velocity. We thus have an increase in resistance 
as the velocity decreases, which has nothing to 
do with the resistance to slip of the solid particles. 

It is similar with the bound water. Here also 
the resistance decreases as the velocity increases. 
In a state of rest the mass held together by 
bound water is essentially elastic, and then 
7=transverse modulus of elasticity X shearing 
stress s. This resistence is overcome when 
motion starts, and at high velocities is entirely 
destroyed. It is not to be assumed that the water 
held in the capillaries acts essentially differently 
from the bound water in the solvation envelopes. 
This causes the difficulty in experimenting that 
the resistance for a given velocity has a different 
value according to whether it is reached by 
acceleration from a lower velocity, or by re- 
tardation of a higher velocity. In the first case 
it is greater than it should be, in the second case, 
smaller; for the destruction of the structure 
requires time, and so does building it up again. 
It is also not to be forgotten that the sludge is 
constantly changing in consequence of biological 
and chemical processes as well as by mechanical 
stirring. 

The flow characteristics will vary according to 
the proportions in which the four components are 
present in the entire mass. 

Corresponding to the definition viscosity 
u=1/fluidity, ¢ for simple viscous liquids, ac- 
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Fic. 2. The five types of flow of plastic bodies. Velocity 
distribution over bore of tubes. v and 7 are velocity and 
gradient in the tube. s and T are the shear and torque 
in the rotation viscometer. The lower part of figure 
shows variation in velocity across the pipe section. 


cording to Bingham we should have for the 
resistance of a plastic material, u’ = 1/mobility, ¢’. 
Since v=Q/F is the average velocity in the 
section of the pipe, ¢’ should be the reciprocal 
value of the average resistance of the material in 
the pipe. (Q=quantity discharged in unit time; 
F=section of pipe.) 

With a shearing stress 7 less than a certain To 
in consequence of the friction or interlocking of 
the solid particles, and also when the jelly 
strength is high enough also on this account, no 
flow will take place or at most an elastic deforma- 
tion. (Regime I of Figs. 2 and 3.) As 7 increases, 
a layer is sheared off in the plane of highest 
stress, that is, near the walls of a pipe, or in the 
case of a rotation viscometer, near the surface of 
the rotating body or of the inner cylinder. The 
mass then moves like a plug through the pipe, 
without further breakdown in structure (regime 
II in Figs. 2 and 3).'* '4 In the rotation viscom- 
eter the inner cylinder slips along a plane of 
shear in the otherwise undisturbed sludge. In 
this case therefore, there is the shearing of one 
layer of liquid between two solid bodies. When 
v increases, the resistance to flow remains con- 
stant for a while (dv/dJ=constant). Finally, 
however, the shearing stress transferred from 
the plane of shear exceeds the resistance of the 
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Fic, 3. Mobility ¢’ =1/u plotted against shearing stress r. 


adjacent layer, and so on; more and more layers 
(in a pipe, concentric rings) come into motion 
relative to the core remaining fixed at the axis 
of the pipe, i.e., a core of decreasing size runs 
more and more ahead of the surrounding layers. 
This core, really a rod, considered as a whole, 
could never, theoretically, lose its solidity en- 
tirely, since r=0 at the axis of the pipe, but it is 
carried along. Next regime III of the flow curve 
shows a gradual increase in mobility and the 
curve must approach asymptotically to an 
inclined straight line; as an actual fact it coin- 
cides with this above a certain point. (Fig. 2, 
regime IV, curve a). The straight line of regime 
IV is applicable practically and theoretically to 
the rotation viscometer, since between two con- 
centric cylinders having relative motion, the 
shearing stress or the velocity gradient never 
becomes zero. The intersection of the straight 
line extrapolated backwards, with the axis of 
abscissas, gives the rigidity (Bingham’s “yield 
value’’) of the material; its dimensions are those 
of a stress, g/cm’. 

If the sludge has no third constituent, no solid 
lyophobe component, which has a certain 
minimum resistance, but also no remaining 
structure in the lyophile parts which is not 
affected by the shearing, if therefore the sludge 
has only an initial jelly strength which can be 
completely destroyed, then in regime IV of the 
flow curve J, it will be represented by a straight 
line which passes through the origin when extra- 
polated backwards, which corresponds to a true 
Poiseuille flow. 

Finally turbulence begins in regime V, within 
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which, with increasing distance from the critical 
velocity, the characteristic differences of the 4 
flow constituents disappear. (Concerning the ve- 
locity distribution in pipes in turbulent flow.)*: * 

Naturally between curves IVa and IVb there 
is the possibility of curves like 1Vc which, without 
any intermediate straight portion, go over im- 
perceptibly into the curve of turbulent flow, 
concave downward. 

Fig. 3 is obtained if, instead of the values v 
and J, the shearing stress and mobility are 
plotted, as for example they are obtained in the 
rotation viscometer. Since the mobility ¢’ is a 
function—to be sure different for each regime 
of flow—of 7 and v is a multiple of y’, the curves 
of Fig. 2 may be regarded as integral curves of 
Fig. 3. The peculiar maximum value of ¢’ to be 
observed in Fig. 3, curve 6, which corresponds 
to the change in curvature in Fig. 2, curve 3, is 
easily explained. Since a pipe not completely 
full (about 90 percent) will carry more water 
than when entirely full, so also a plastic mass, 
which has a still unsheared core around the axis 
of the pipe, may have a lower resistance than 
when the parabolic velocity profile has been 
formed across the entire pipe section. For with 
the same rate of flow, in the first case the surfaces 
of the layers slipping past one another are 
essentially smaller, while the velocity gradient 
is only slightly larger than in the second case. 

The curves, sketched above, of Figs. 2 and 3, 
are to a certain extent ideal. Only very seldom 
is it possible to distinguish between all five 
regimes and recognize them. With a not homo- 
geneous material, such as sewage sludge, they will 
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Fic. 4. V plotted against pressure P (after Ostwald). 
A is for mercury sulpho salicylic acid sol. B is for a 1 
percent gelatine sol. 





WILHELM MERKEL 


generally blend more or less. If turbulence 
occurs early, regime IV may be left out com- 
pletely, with increasing water content regime II 
becomes III, and also regime I may completely 
disappear (Fig. 4). 

Until now it has been tacitly assumed that 
only such a sludge would be considered as con- 
tained enough water so that pumping it through 
pipe lines would come under consideration. With 
an air-dry, fermented sludge, the constituent 
No. 1 (free water) and No. 2 (slimy material) 
are sO unimportant in comparison with No. 3 
(lyophobe particles) that, in consequence of the 
large friction of these solid particles, flow under 
pressure through the pipes is no longer possible. 
For when lubrication is lacking the friction of 
solid bodies increases with increased pressure 
(when pressed together the surfaces of contact 
are increased by elastic yielding) until finally 
pressure arches are formed in the pipe. These 
arches are first solidified by further increase in 
pressure, and are destroyed only at much higher 
pressures. Not until the sludge is completely 
dried (cured) could conveying through pipes be 
again considered, but then in powdered form 
transported by air, i.e. with air as the vehicle. 

From what has been said there is no prospect 
of setting up a single general equation, to cover 
all the variations and the separate regimes of 
flow of sewage sludge or of plastic materials in 
general. It is therefore necessary to start with 
data tables giving supplements for the most 
varied effects, such as are available for draining 
land, in the field of cultivation. A diagram ac- 
cording to the available experimental data is 
shown in Fig. 5. There is still need of the coop- 
eration of colleagues to determine the coefficients 
more exactly and to carry the work further until 
it suffices for the manifold relations of practice. 

Now if no general equation of flow can be 
obtained, there may yet be found a suitable 
formula for a restricted field. Reiner? derived 
such a one in seeking a general law. In the 
regimes of flow III and IV, which with different 
kinds of sludge generally occur between 80 and 
88 percent water content, it may be used with 
good approximation to convert to other ¢, v 
and J. It is derived from the relation r=0 


+ ydv/dr and reads, 
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Fic. 5. Pressure loss of sewage sludge in pipes 20 cm in diameter. 
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Fic. 6. Sewage sludge from a leather factory at 20°C 
(viscometer tests). Mobility ¢’ and stiffness yp’ (physical 
units) plotted against falling weight G in grams. The two 
lowest graphs are marked ‘‘unworked” and “‘worked,”’ i.e., 
unstirred or stirred before testing. 
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Here ¢ is to be taken in cm and yz in the technical 
unit g-sec./cm’; # has the unit g/cm? and y is 
g/cm*. Without the quantity in brackets we 
have the Poiseuille equation, which is thus a 
special case of this general equation. The third 
term in the brackets can be neglected on account 
of its small value above a limiting “structure.” 
The curve therefore goes over to a straight line. 
In the curves calculated from the above equation, 
their asymptotes are parallel to the Poiseuille 
lines with the same yu. In the regime where the 
third term can be neglected, the equation may be 
written more simply, 
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of Eq. (2), as compared with the Poiseuille 
straight lines, v=?r?>Jy/800u, becomes f=J,-5 
= 2.6673/ry. Two known points suffice for con- 
version in the IVth regime of flow. 

With increasing fermentation, naturally with 
the same water content, » decreases and # 
increases. (See Fig. 5 with the steeper curves for 
the completely fermented summer sludge as 
compared with the flatter ones for the only half 
fermented sludge of low winter temperatures.) 

For comparison with the experimental results 
with sewage sludges reference will be made to 
publications which treat of other plastic materi- 
als and show the same results. Scott Blair'* 
shows the flow curve of a loam paste which 
shows very beautifully the above mentioned 
regimes of flow. For loam, as a mixture of lyophile 
clay and lyophobe sand, is very similar to sewage 
sludge. Further the curves for zinc sulphide in 
a-bromonaphthalene! show an agreement with 
Fig. 5. 

As an especially instructive example of the 
results obtained with the Kaempf viscometer, 
curves are shown in Fig. 6, which were obtained 
in tests of sewage sludge from a leather factory. 
Noteworthy is the evident occurrence of the II 
regime of flow (see Fig. 3). Since the thixotropic 
stiffening of the leather sludge was not great and 
the stiffening time was relatively long, the 
separate values show only slight deviations on 
repetition. Further details are given in the 
paper of the author above referred to. 

In addition to the direct significance which the 
flow curves obtained have for the conveyance of 
sewage sludge, advantages are to be expected 
from applying the theoretical results to other 
plastic materials. In all industries where plastic 
materials are prepared, mixed, pumped and 
worked, flow curves should be prepared for the 
material in question, and it will surely be seen in 
many cases that the efficiency may be consider- 
ably increased by changing the velocity of mixing 
or of flow. Especially in the IV regime of flow a 
certain increase in driving force corresponds to 
quite a considerable increase in output. If a 
material is thixotropic, the increase in efficiency 
by increase in motor capacity is even more 
striking. 

In this connection may be remembered the use 
of grout, the preparation of road foundations 
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with tar or bitumen, pumping paper stock and 
molasses, and the wet transportation of any kind 
of refuse. 

In summary it may be said, if any kind of 
thick material is to be investigated as regards 
its flow characteristics, a clear distinction must 
be made between different materials: 

1. Clean or slightly soiled liquids. It is first 
to be determined whether or not the flow is 
turbulent at the expected velocity. In turbulent 
flow the well known formulas for water are to 
be applied with suitably changed coefficients; 
the roughness of the walls is to be taken into 
consideration. 

If laminar flow is to be expected, Poiseuille’s 
law applies. Only the viscosity has to be deter- 
mined; a moderate roughness of the walls plays 
no part. 

With increasing contamination of the liquid, 
the validity of the above methods of calculation 
disappears and it is to be determined whether 
and to what extent the suspended solid particles 
attract and bind the liquid to them (lyophile) or 
not (lyophobe). 

2. With increasing content of lyophile par- 
ticles, the viscosity uw increases very rapidly and 
thixotropic phenomena are also to be expected. 
The critical velocity increases, below it there may 
be a true laminar regime (a straight velocity- 
pressure drop line through the origin). At lower 
velocities, greater resistance, and under certain 
circumstances, a yield value may appear. 

3. Increasing addition of lyophobe material 
does not necessarily have a large effect on the 


viscosity, so long as it does not become so great 
that capillary action occurs between the closely 
packed solid particles. Here appears a flow re- 
sistance &. Higher concentrations all have a 
yield value below which no flow is possible. At 
high velocities the inner friction can never 
entirely disappear. There is only over a certain 
range a pseudolaminar flow, i.e., a flow-pressure 
drop line which does not pass through the 
origin. If the amount of the solid lyophobe par- 
ticles is greater than the amount of liquid, the 
friction of solid bodies comes to the foreground. 
Turbulence is no longer possible and friction 
depends more on the pressure than on the 
velocity.‘ 

The author carried out this investigation at 
the suggestion and with the support of Professor 
E. Neumann, Techn. Hochschule, Stuttgart. 
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